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ABSTRACT: We resurrect a standard construction of analytical mechanics 
dating from the last century. The technique allows one to pass from any dynam- 
ical system whose first order evolution equations are known, and whose bracket 
algebra is not degenerate, to a system of canonical variables and a non-zero Hamil- 
tonian that generates their evolution. We advocate using this method to infer a 
canonical formalism, as a prelude to quantization, for systems in which the naive 
Hamiltonian is constrained to vanish. The construction agrees with the usual re- 
sults for gauge theories and can be applied as well to gravity, even when the spatial 
manifold is closed. As an example, we construct such a reduced Hamiltonian in 
perturbation theory around a flat background on the manifold x R. The result- 
ing Hamiltonian is positive semidefinite and agrees with the A.D.M. energy in the 
limit that deviations from flat space remain localized as the toroidal radii become 
infinite. We also obtain closed form expressions for the reduced Hamiltonians of 
two minisuperspace truncations. Although our results are classical they can be 
formally quantized to give the naive functional formalism. This is not only an 
effective starting point for calculations, it also seems to provide a formulation of 
the quantum theory which is non-perturbative, at least in principle. The marriage 
we advocate between the old technique and canonical quantization seems to have 
profound implications for quantum gravity, especially as regards the conservation 
of energy, statistical mechanics, and the problem of time. 
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1. Introduction 

It is traditional to develop the canonical formalism from a Lagrangian. One first 
identifies variables whose Poisson bracket algebra is canonical, then the Hamiltonian is 
constructed and used to generate first order evolution equations. However, the imposition 
of constraints can result in a system of dynamical variables whose bracket algebra is not 
canonical and whose time evolution, while completely determined, is not generated by a 
known Hamiltonian. This possibility is especially relevant to systems in which gravity 
is dynamical and the spatial manifold is closed because then the naive Hamiltonian is 
constrained to vanish [1] . The ubiquity of canonical quantization and its evident failure in 
this case has led to much puzzlement over how these perfectly acceptable classical systems 
should be quantized [2]. We wish here to advocate a straightforward solution: impose 
the constraints to obtain directly a set of first order evolution equations and a bracket 
algebra for the reduced variables, then construct a canonical formalism which reproduces 
this system. The resulting reduced canonical formalism can then be quantized as usual. 

Our motive in suggesting this step is the highly unsatisfactory situation which currently 
prevails in the formulation of canonical quantum gravity. Those who approach this field 
from the perspective of other disciplines have long been frustrated by its failure to obey 
certain obvious correspondence limits. There are at least four problems of this nature: 

(1) The Paradox of Second Coordinatization — In classical gravitation we know that fixing 
the lapse and shift determines the coordinate system up to deformations of the initial 
value surface. We also think we understand pretty well how to infer physics from 
the behavior of the metric field in these coordinates. In quantum gravity on a closed 
spatial manifold we are told that even after fixing the lapse and the shift we must 
determine what time is all over again. We are also told that physical information can 
only be gleaned from observables which are manifestly coordinate invariant, that we 
don't possess any such observables, and that even if we did it would not be possible to 
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compute their expectation values because an inner product cannot be given until after 
the meaning of time has been clarified. How can changing % from zero so thoroughly 
confuse the way we infer physics from the field variables? 

(2) The Paradox of Dynamics — A similarly striking issue concerns the limit in which 
gravity becomes non-dynamical. Since the Hamiltonians of pure matter theories are 
not typically zero, even on closed spatial manifolds, these theories can be quantized 
canonically. But when gravity is made dynamical on a spatially closed manifold it is 
asserted that we no longer know even what the inner product is, no matter how weak 
the gravitational interaction might be relative to other forces. Note that the alleged 
difficulty is unrelated to the non-renormalizability of Einstein's theory; it would occur 
even in an ultraviolet finite theory of quantum gravity. How can changing Newton's 
constant from zero affect the basic structure of quantum mechanics? 

(3) The Paradox of Topology — Although topology is not a continuous parameter the 
comparison between manifolds with differing topologies is disquieting. For spatially 
open manifolds there is no obstacle to canonical quantization because the naive Hamil- 
tonian can be non-zero [3,4]. Of course there remains the highly non-trivial problem 
of finding a theory whose dynamics are consistent, but there is no confusion about 
basic issues such as the meaning of time or how to define inner products. When the 
spatial manifold is closed it is asserted that we no longer understand these issues. Yet 
there is no local experiment which can distinguish a closed space from an open one, 
as witness the fact that either possibility might describe our own universe. How can 
removing a point from the other end of the universe affect our ability to provide a 
quantum mechanical description of local observations which have a vanishingly small 
probability for even being in causal contact with the boundary? 

(4) The Paradox of Stability — Finally, there is the peculiarity encountered in passing 
from pure quantum mechanics to statistical mechanics on a spatially closed manifold 
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when gravity is dynamical. Since the energy is constrained to be zero it follows that 
all states are degenerate. For example, the total energy needed to add a particle-anti 
particle pair to any state is zero because the negative gravitational interaction energy 
cancels the positive energy of the pair. Simple considerations of entropy seem then to 
require that the microcanonical ensemble should be concentrated around what would 
otherwise be thought of as very highly excited states. To see this note that for every 
"empty" state there are a countably infinite number of degenerate states which contain 
a single particle-anti particle pair moving apart with various momenta. There are even 
more states which contain two pairs, and more yet which contain three, etc. So what 
prevents such a universe from evaporating into a maelstrom of pairs? 

We resolve the paradox of second coordinatization by denying that there is any fun- 
damental distinction between classical and quantum measurement beyond that usually 
imposed by the uncertainty principle. What time means in quantum gravity is determined 
by fixing the lapse and the shift, just as in the classical theory. Any question that can 
be answered classically by studying a functional of the metric with fixed initial value data 
can be studied as well in the quantum theory using the expectation value of this same 
functional in the presence of a state whose probability is concentrated around the classical 
initial value data. We can base observations on a non-invariant such as the metric because 
the gauge has been fixed: any quantity becomes invariant when it is defined in a unique 
coordinate system. 

We resolve the paradoxes of dynamics and of topology by showing that for any con- 
strained system which can be reduced — that is, for which the residual gauge freedom can 
be fixed and the constraints imposed — there is a canonical choice of reduced dynamical 
variables and a non-zero Hamiltonian which generates their time evolution. In fact there 
are many such Hamiltonians, each corresponding to a different identification of canonical 
variables in the reduced dynamical system. For gravity these Hamiltonians seem to have 

4 



the property that the closed space versions go over to the known open space ones when 
the coordinate volume is taken to infinity in such a way that the initial value data are only 
locally disturbed from a background which obeys the appropriate asymptotic conditions.* 
When gravity is coupled to matter the non-zero reduced Hamiltonians seem to agree with 
the corresponding pure matter Hamiltonians in the limit that Newton's constant vanishes. 

Although the Hamiltonians we are discussing generate time evolution they are not 
generally conserved. Outside of a few special cases the generally conserved energy really 
is zero, so the universe can potentially suffer from the instability described in (4). We see 
a mostly empty universe today because the balance implicit in the H = constraint is 
maintained by means of global, negative energy modes, and the rate at which these modes 
can be excited is highly suppressed by causality and by the weakness of the gravitational 
interaction. There has simply not been enough time for the universe to evaporate into 
pairs. Depending on the topology and the causal structure of the initial state there may 
never be enough time for this to happen. 

Section 2 presents the construction by means of which any known system of first order 

evolution equations and (not necessarily canonical) brackets can be used to identify a set 

of canonical variables and a Hamiltonian that generates their time evolution. This is not 

original work; the construction was given about a century ago [5,6]. What does seem to be 

original is the idea of applying it generally to gravity. In taking this step we are following 

the work in 2 + 1 dimensions of Moncrief [7], Hosoya and Nakao [8], and Carlip [9]. We 

believe that their constructions are special cases of the general technique; we believe the 

same is true of the lovely formalism for open, asymptotically flat spaces which was devised 

by Arnowitt, Deser and Misner (A.D.M.) [3, and references therein]. 

* For [3] this background would be flat space, for [4] it would be the de Sitter geometry. 
Note that topological obstructions — for example, the inability to impose a flat metric on 
— need have no physical relevance because we might be able to use half of the infinite 
coordinate volume to approach the asymptotic geometry of the open manifold and the 
other half to deviate so as to reconcile the topology of the closed manifold. 
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In section 3 we apply the method to a standard gauge in scalar electrodynamics. 
This is difficult to do in gravity because one can seldom obtain explicit solutions for the 
constraints. However, the method is simple enough to carry out in perturbation theory or 
when the dynamics is suitably truncated. In section 4 we give a perturbative construction 
for a flat space background on the manifold x R; in section 5 we apply the method to 
several models of minisuperspace. We show in section 6 that if the unreduced theory exists 
then quantization along the lines we advocate results in the usual functional formalism. 
In this form the technique is not only simple to exploit for the purposes of perturbative 
calculations, it also provides a definition of the theory that is valid beyond perturbation 
theory, at least in principle. Section 7 summarizes our resolutions to the aforementioned 
paradoxes and discusses the implications our view has for statistical mechanics and the 
conservation of energy. 

2. The Construction 

Consider a set of 2N dynamical variables, {v l {t) : R — > R 2N }.* Suppose that the 
time evolution of these variables is determined by first order equations of the form: 

v\t) = f%t) (2.1) 

where the / z, s are known functions of the v l, s and possibly also of time. We do not assume 
that the v l 's are necessarily canonical but rather that they obey the following bracket 
algebra: 

{v i (t)y(t)} = J i 3(v,t) (2.2) 

where the bracket matrix 3 l] is antisymmetric, invertible and obeys the Jacobi identity. 
We shall denote the inverse bracket matrix by the symbol, "Jjj": 

3^3^ = 8^ = 3^ J 'ji (2.3) 

To emphasize that the dynamical variables are functions of both the canonical positions 
and their conjugate momenta we have adopted the neutral symbol, 'V," for "variable." 
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Taking the time derivative of (2.2) and substituting (2.1) reveals the following relation: 

J ik f i , k -J ik f i ,k = J ij . k f k + d -^- (2-4«) 

where a comma denotes differentiation. The analogous relation for the inverse is: 

d J ■ 

J ik f ,j ~ Jjkf ,i = ~ J ij,k f -Qj~ ( 2Ab ) 

Since brackets obey the Jacobi identity we have the following differential relation for 
the elements of the bracket matrix: 

J U J jk £ + P l J k \ ( + J M J ij £ = (2.5a) 

The analogous result for the inverse is: 

Jij,k + J jk,i + J ki,j = ( 2 - 5& ) 

This last relation is reminiscent of the Bianchi identities in electromagnetism and allows 
us to write the inverse bracket matrix as the curl of a "vector potential," (f)f 

Jij = 4> 3 ,i ~ 4>i,j (2-6) 

Just as in electrodynamics fa is undetermined up to the gradient of a scalar function; 
similarly, one may need several coordinate patches if the dynamical variables map into 
some space other than R^ N . A convenient representation for the vector potential is: 

( j) i {v^t) = -l drr Jij(vQ + rAv,t) Av j (2.7) 
J 

where Vq is any fixed point and Av l = v l — Vq. Note that under a change in the dynamical 
variables J lJ transforms contravariantly while Jij transforms covariantly; fa transforms 
covariantly up to the addition of a gradient. 

The simplest possibility for reconstructing a Hamiltonian is that the evolution equa- 
tions are integrable. The condition for this is that the bracket matrix should be free of 
explicit time dependence: 

dJ lJ 

3H(v, t) 3 v(t) = {v\ H} iff — — = (2.8) 
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For the proof note that, if it exists, the Hamiltonian is determined up to a function of time 
by the equations: 

Since we must have Hjj = Hji the integrability condition is: 

= (J ik f k ) . ~ (J jk /*) . (2-lOa) 
= " ( J ki,j + Jjk,i) f k + (Jik f k j - f k ,i Jjk) (2-106) 

where we have used (2.46) and (2.56) to reach the final line. The theorem is proved upon 
noting that the bracket matrix is invertible. When the integrability condition is met it is 
simple to check that: 

H(v,t)-H(v ,t) = drAv'Jijivo + TAv^ffvo + rAv,^ (2.11) 
J 

gives an explicit representation for the Hamiltonian. 

To identify canonical variables in the local neighborhood of any point Vq we construct 
2N functions Q a (v , t) and P},(v, t) (a, 6 = 1, 2, ... , N) which are instantaneously invertible 
at any time and which obey: 

dO a 

Mv,t) = Pa(v,t)-^ r (v,t) (2.12) 



Note that invertibility implies: 

6 b \ _ I W dv* m \ (2.13a) 



6 b J dfadvL dPadv 

a 7 v dv i dQ b dv* ~dP b 



dv l 8Q a , dv % dP a 



The problem of showing that such functions exist is known as Pfaff's problem in honor 
of the German mathematician J. F. Pfaff [10]. This problem was solved long ago by G. 
Frobenius and J. G. Darboux [11,12,13]. That the transformation: 

q a (t) = Q a (v(t),t) , p a (t) = Pa(v(t),t) (2.14) 
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gives canonical variables is straightforward to see. First, substitute (2.12) into (2.6) to 
obtain: 

= dPadQ^_dQ^dPa 

lj cV dvi cV dvi [ ' ' 
Now use this relation and (2.13a) to show that the bracket matrix "raises" indices on the 
transformation matrices: 

dv i ■ ■ dv k 

w= J3J *w (2 - 16a) 

Jl3 ( dP^d^_d^dP^dv^_ 

\dvi dv k dvi dv k ' dQ a 1 ' ' 

= jijd iB (2 - i6c) 

% •/"./ / ,^ (2.17a) 

8P a ]k dP a V ' 

_ fdP^d^__d^dP^\d^_ 

\dv3 dv k dvi dv k ) dP a { ' ' 

= ~J ij % (2.17c) 
Finally, take the brackets using (2.2), (2.13a), (2.14), (2.16) and (2.17): 

{8 .,rt = ^J«^ = -^^=C (2.18a) 

XH ,H s dv l cV dv l dP b v ' 

{ S, K) = WjVm = °r°*=t} (2.186) 

iy Fbf cV dvi dv % dQ b a v ' 

{Mh) = ^r^ = ^^-« (2 . 18 c) 

A constructive — and slightly different — procedure for finding canonical coordinates is 
credited to the French mathematician Darboux and can be found in the text by Arnold 
[14]. The older derivation is discussed in the text by Whittaker [15]. 

It is possible to give explicit formulae for N = 1. For this case we can write down the 
inverse bracket matrix: 

t t) = (i T i < 219 » 
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and we further know that the function J 12 (v, t) is non-zero. The most convenient "gauge" 
for the vector potential is an axial one: 

MvW,t) = j Q 1TW ^ J) (2.20a) 

<fo(v 1 ,v 2 ,t) = (2.206) 

Comparison with (2.12) reveals that Q cannot depend upon v 2 and that the choice of Q 
uniquely determines P. The simplest choice is: 

Q(v 1 ,v 2 ,t) = v 1 (2.21a) 

Note that although the general construction is only local — and so must be built up in 
patches — the formulae for N = 1 are valid globally if the domain of (v 1 , v 2 ) is R 2 . 

Of course the mapping to canonical variables can only be unique up to canonical trans- 
formations — even in the Lagrangian formalism. This shows up slightly in the integration 
constants which come from solving (2.12); for the case of N = 1 above we had the freedom 
to choose Q to be any function of v 2 and t which was instantaneously invertible on v 2 . By 
far the larger source of variation is the freedom to add the gradient of a scalar to ^. From 
the instantaneous relation: 

fadv* = Pa dq a (2.22) 

we see that the "gauge transformation," <^>j — > <pi—di9, induces a canonical transformation 
whose generating function is —9. 

Once the brackets are canonical they are also independent of time, so we can get a 
Hamiltonian in the new coordinate system by using (2.11). The evolution equations in 
these coordinates are: 

q a (t)=f a (q(t),p(t),t) (2.23a) 
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Pa{t)=f a (q{t),p{t),t) (2.236) 

where we define: 

/« (Q(v, t), P{v, f ), f) = ^(v, t) f(v, t) + t) (2.24a) 

/ \ 8P ■ BP 

f a [Q(v, t),P(v, f), f) = -^(v, t) f(v, t) + ^(v, t) (2.246) 

Substituting into (2.11) gives: 
H(q,p,t) - H(q ,p Q ,t) 

= J dr ^Ap a f a (q + TAq,p + TAp,t) - Aq a f a (q + T Aq,p Q + rAp,t^ 
where (qo, Po) is any point in phase space and we define Aq a = q a — qfi and Ap a = p a — poa- 



Note that the Hamiltonian (2.25) only generates the evolution (2.23) of the canonical 
coordinates. Its bracket with v l (t) does not generally give v l (t). If we call the inverse 
transformation V l (q,p, t) then the relation: 

v i (t) = V i (q(t),p(t),t) (2.26) 

implies that the original dynamical variables acquire only a portion of their time depen- 
dence from that of the canonical coordinates: 

dV i dV i dV i 



r ■ i dV l 

{v t (t),H} + — (2.276) 



We have already seen that the final term in (2.276) must be non-zero if the J tJ, s contain 
explicit time dependence. 

Far from being a problem, this extra source of time dependence is a blessing. It is what 
keeps the evolution of the original dynamical variables fixed if we alter what is meant by 
the canonical variables. As long as there is a clear procedure for inferring physics from 
the v l (tys it does not matter that we can choose the canonical variables differently. This 
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is a subtle point and a crucially important one. It is also counterintuitive because any 
choice of canonical coordinates does offer a complete description of the degrees of freedom 
available to the system. The point is that what these degrees of freedom mean physically 
changes as we change our representation of them. 

The possibility for this sort of subtlety is always present in physics because time de- 
pendent canonical transformations can be made in any classical theory and, with proper 
attention to operator ordering, on the quantum level as well. For example, in perturba- 
tive quantum field theory there is a well known transformation, called as the "interaction 
representation," which takes one to a set of fields whose time evolution and commuta- 
tion relations are free. Consider the theory of a scalar field, (f)(t,x), whose action has the 
following form: 

S[<f>} = J d 4 x{^cf>d^- \m 2 ^ 2 } + SM (2.28) 

where Sj[4>] is an ultralocal interaction. Since Sj[4>] is free of derivatives the momentum 
canonically conjugate to <j)(t, x) is just its time derivative. The only non-vanishing equal 
time bracket is: 

{<P(t,x)^(t,y)} = 5 3 (x-y) (2.29) 

The field equation: 

(□-m 2 )^ + ^M = (2.30) 
can be integrated to give a new field which obeys the Klein-Gordon equation: 

□ — m z d(p 

If we define the inverse of □ — m 2 so that it and its first derivative vanish at t = then we 
have also: 

$(0,x) = 0(0, f) (2.32a) 

6(0, f) = 0(0, f) (2.326) 
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It follows that <&(£, x) not only evolves like a free field, it also obeys the same bracket 
algebra. We do not conclude that all such theories are free, even perturbatively, because 
we insist upon inferring physics from the original variables. The only difference between 
this standard situation and the previous discussion of this section is that even the original 
scalar is canonical whereas we do not make this assumption about the v l, s. 

Why the imposition of constraints tends to result in a non-canonical bracket matrix is 
obvious to anyone who has ever constructed a Dirac bracket [16] . Even though the gauge 
and constrained variables are irrelevant to physics, the unreduced canonical formalism is 
an artificial construct in which the unphysical variables play an essential role in keeping 
the brackets canonical under time evolution. After reduction these variables are no longer 
independent and this changes the bracket algebra. Unless what we call the reduced vari- 
ables are chosen very carefully they will not have a canonical bracket algebra. This does 
not mean that other choices are "wrong." In the sense of providing a complete description 
of physics a non-canonical set of variables can be as "right" as a canonical set. Non- 
canonical reduced variables can even be highly preferred on account of bearing a simpler 
relation than any canonical variables to the original degrees of freedom from which physics 
is inferred. 

In the next section we use scalar electrodynamics to provide an explicit example of how 
reduction can yield a non-canonical bracket matrix. For now we can understand the issues 
in a simple but somewhat contrived fashion through a model of two coupled oscillators 
whose Lagrangian is: 



The resulting canonical formalism will serve as our model for an unreduced canonical 
formalism, in spite of the fact that the Lagrangian possesses no continuous symmetries 
and gives rise to no constraint equations. We shall make up for the absence of constraints 
by imposing them ad hoc. 




(2.33) 
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It is straightforward to show that the canonical variables of the coupled oscillator 
system have the following time evolution: 

Ql(t) = \ + %) cos (-art) + \(q\ - q 2 ) cos(Jart) 

2V [ 3 ! ! (2.34a) 



Pl(t) = --^ruo(qi + (f 2 ) sin (-art) - ^mu(qi - qo}j sin(^art 

+ ^ (pi + P2j cos (2^) + 2 (pi _ cos (2^ 

92 (*) = J(<fl +$2) cos(^wt) + -(-$1 + 92) cos(Jart) 

2V 13 1 1 ( 2 - 34c ) 



P2(t) = -\mu(qi + q 2 ) sin (-cut) - Jmwf-gi + q 2 ) sin (-art) 

4 V x y 3 1 1 (2 - 34d) 

+ 2(^1 +^2) cos (-a;*) + - (-pi + p 2 ) cos(-u;^ 



In these formulae we have denoted the initial values by a hat. It is these initial values 
that are the independent degrees of freedom of this system; in particular, only they have 
independent Poisson brackets. It is straightforward to verify that if the initial variables 
have Poisson brackets then the time evolution given in (2.34) preserves this feature. That 
is, the only non-zero equal time brackets are: 

= 1 = {q 2 (t),P2(t)} (2.35) 

We emphasize that the bracket algebras for different times are not independently specifi- 
able. Since the variables at time t are uniquely determined in terms of their initial values, 
the Poisson bracket algebra at time t is determined by time evolution and by the initial 
bracket algebra. 

The model at this stage should be thought of as analogous to temporal gauge QED 
or Yang-Mills, or to synchronous gauge gravity. A constraint in these models is a relation 
between the initial value variables, and such constraints are reduced by imposing a gauge 
condition on the initial value surface. It is only in exceptional cases that we can find 
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gauge conditions which are preserved under time evolution. Initial value gauge conditions 
always imply some relation between the later canonical variables, but very seldom the same 
relation. Unless we choose what we call the variables of the reduced theory to compensate 
for this change, the reduced brackets will become non-canonical because at any instant 
they are the Dirac brackets associated with different gauge conditions. 

Let us suppose that the reduction of our oscillator model is accomplished by setting 
§2 = = P2- Let us further suppose that we take as our reduced variables, v^(t) = q\(t) 
and v 2 (t) = pi(t). Since the constraint and its conjugate gauge condition affect only the 
initial values we see from (2.34) that the reduced variables have the following evolution: 

v l (t) = -oi cosf-cut^ H — oi cos( -ut) H m sinf-cut^ H pi sin( -ut] (2.36a) 

v 1 2 yi V2 / 2 yi V2 / 3miu V2 / m/ 1 V2 / v ' 

2/ x s _ /3 \ i ^ /i \ i_ /3 \ i^ n \ , 

v (t) = —-mu)q\ sm[-ujt\—-mujqi sin I -u>t I + -p\ cost -cot ) + -pi cos I -cut) (2.366) 

4r \ Zt ' 4r \ Zj / Zj \ Z ' Z \ Z ' 

It follows that the only non-zero equal time Poisson bracket is: 



cos -ut 



- COS" I ... / 



(2.37) 



Although it is canonical at t = 0, it is not so later on, and it can even pass through zero! 

If we ignore the fact that the bracket matrix becomes non-invertible at cot m (2n + 
l)n ± .127T then the construction gives the following canonical variables: 



P(t) = 



q(t) = v\t) 



4 cos 2 (jz^tj — cos 2 (utj 



The associated Hamiltonian is: 



H(t,q,p) = ^cos 2 (^iut) - icos 2 (W) 



1 2 2 
+ -muj q 



V 

2m 

cos J ( tjoot^j + -j cos- [ut 
cos 2 ^ryo;tj — |cos 2 ^o;tj 



(2.38a) 
(2.386) 



(2.39) 
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Note that although this Hamiltonian generates time evolution for q(t) and p{t), it does 
not do so for either the reduced variables or for the original ones. The variable v 2 (t) 
acquires additional time dependence from relation (2.386), and the original variables have 
the following expressions in terms of q(t) and p(t): 



Ql(t)=q(t) 



Pl(t) = 



4 ofl \ 1 o, 

- cos ( -ut I cos I ut 

3 V 2 / 3 



p(t) 



Q2(t) 



P2(t) = 



-2sin(£urt) sin(furt) 4. 2 /l \ / \ p(t) 
7^ — ^ v 7 \ q(t) + -sm 2 [-ut) sin (ut) — 

lra C 2 i,,rf _,™2 ,rf 3 \2 J \ / mu 



4 cos^ ( r^ut ) — cos z ( ut 



— 3cos 2 (^o;t ) sin [art 



mwg(f) + | sin (-art) sin (-art) p(t) 
4 cos 2 Uart ) -cos 2 (art ) 3 V2 / V2 / 



(2.40a) 
(2.406) 

(2.40c) 
(2.40d) 



Although the preceding canonical formalism is valid, it is not simple because the re- 
duced variables were badly chosen. A much more convenient choice is: 

v 1 (t) = qi (t) + q 2 (t) = qi cos (-ut) + — pi sin (-a;*) (2.41a) 

v 2 (t) =Pi(t) +P2(t) = --muq\ sin(^art) +pi cos(^art) (2.416) 
The resulting reduced bracket algebra is canonical: 

{v 1 (t),v\t)} = l 

and the Hamiltonian is time independent: 



(2.42a) 



H = — (v 2 ) 2 + -mu 2 [v l 
2m V J 8 



(2.426) 



This Hamiltonian generates the evolution of f 1 (t) and v 2 (t) but not that of the original 
variables. They acquire additional time dependence through the relations: 



qi (t) = 2 cos 2 (^art) — cos 2 (art) ^(t) — ^sin 2 (^ort) sin (cut) 



v 2 (t) 



mu 



(2.43a) 
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Pl(t) = cos 2 (^-ut^j sin^cutj mujv^it) + 



r2 o(\ \ 1 o, 
- cos ( -cut ) + - cos ( cut 



t; 2 (t) 



©(f) 



■ 2/1 
2 sin ( -u)t 



sin 



^cutj t; "*"(£) + — sin 2 ) 



mcu 



P2(f) = — cos 2 ^o;tj sin^cutj muv 2 {t) + ^sin 2 ^u;tj + ^sin 2 ^u;tj f 1 (f) 



(2.436) 

(2.43c) 
(2.43d) 



Another convenient choice is: 



(2.44a) 
(2.446) 



(2.45a) 



vl (t) = Ql(t) - = qi cos(^cut) + — pi smf^cvt 

f 2 (f) =Pl(*) ~P2(t) = -^muqi sin^cut) +pi cos^-cu^ 
As before, the resulting reduced bracket algebra is canonical: 

{v 1 (t),v 2 (t)} = l 

and the Hamiltonian is time independent: 

This Hamiltonian generates the evolution of v 1 (t) and v 2 (t). The original variables acquire 
additional time dependence through the relations: 



(2.456) 



11 {t) = ^ cos2 {\ ut ) Jr \ cos2 { ut ) ^(t) + ^ sin 2 (^ut^j sin(ut 



v 2 (t) 



P\(t) = — cos 2 ^o;tj sin(ut^ muv l (t) + 2cos 2 ^u;tj — cos 2 ^u;tj v 2 (t) 



92(f) = ~~ g sin 2 (2^) ~~ 3 s ^ n2 ( a; *) 1,1 3 s ^ n2 (2^) s ^ n ( a; ^ 



^ 2 (f) 



mcu 



(2.46a) 
(2.466) 
(2.46c) 



p 2 (t) = -cos 2 (^cut) sin(W) muv l {t) + -2 sin 2 + sin 2 (W) w 2 (t) (2.46a 1 ) 

Three points deserve mention. First, although the latter two formulations are signifi- 
cantly simpler than the first, all three give precisely the same evolution for qi(t) and Pi(t). 
As long as we infer physics from these variables the seemingly different theories are in fact 
identical. Physicists are so conditioned to regard the Hamiltonian as an observable that 
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this point can not be overemphasized. A glance at the first Hamiltonian (2.39) suggests a 
harmonic oscillator with time dependent mass and frequency: 

TTl 

mit) = 4 o/l \ 1 of ^ {2A7a) 
S COS / ikUlt — 4 cos z ut 



u 2 (t)= cos 2 (7^) + J cos 2 (ut) u 2 (2.476) 



The second Hamiltonian (2.426) is that of an oscillator of mass m and frequency |u;, while 
the third Hamiltonian (2.456) is that of an oscillator with mass m and frequency In 
fact the system is none of these things, it is rather a set of coupled oscillators which has 
been subjected to a constraint. This only becomes apparent by ignoring the Hamiltonian 
and concentrating instead upon the ft(t)'s and the pj(£)'s. 

Note that although the reduced canonical variables provide a complete and minimal 
representation of the system's dynamical degrees of freedom, there is nothing wrong with 
using as observables the overcomplete representation provided by the original variables. 
These will include some pure gauge degrees of freedom and some constrained ones. Of 
course the former are unphysical but the latter encode perfectly valid information, even if 
this information can be recovered from the reduced variables. The longitudinal electric field 
offers a familiar example. We really do not need this quantity since it can be recovered 
from a knowledge of the positions of all the charges, but there is nothing wrong with 
regarding it as an observable since it can be measured using the Lorentz force law. 

The second point is that the reduced canonical variables do have a role to play. In the 
quantum theory they tell us how to describe states, and how the original variables act as 
operators upon these states. For example, suppose we quantize the first representation of 
the coupled oscillator system. A state in the Schrodinger picture position representation 
is described by a square integrable wavefunction, ip(q, t). The time evolution of this wave- 
function is generated by the Hamiltonian (2.39). The operators qi{t) and pi(t) act on such 
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a state through relations (2.40), where q(t)ip(q,t) = qtp(q,t) and p(t)ip(q,t) = —i ^jjf^ . 
Note that even in the Schrodinger picture the observables qi(t) and Pi(t) have time depen- 
dence. In the Heisenberg picture the state is time independent and the evolution of the 
reduced canonical variables, q(t) and p(t), is generated by (2.39). 

The final point concerns energy. We have stated that none of the Hamiltonians - 
(2.39), (2.426) and (2.456) — plays the usual role of the energy as an observable, so what 
does? The answer is the Hamiltonian of the original theory, after the imposition of our 
"gauge condition" and "constraint," <?2 = = P2- 



P2 , 5 2-2 
2¥ + 8 m " gi 



(2.48a) 
(2.486) 



<72=0=P2 

This energy does not generate time evolution for the reduced canonical variables in any of 
the three representations given above. For example, in the last representation we can use 
relations (2.34) and (2.46) to express the energy as: 



E(v\v 2 ,^j = [^ + ^ cos M)j \moJ 2 (V) 



15 



sin(c<;t) u v 2 + 



17 15 1 1 

¥ ~ ¥ cos ^ 



(2.49) 



2m 



From (2.40a) it follows that: 



7-, 1 / \ 1 15 . / s i / s r 17 15 , . 

E, v (t) \ = — — sm(ut)uv (t) + — — cos(utf) 

Jo L 8 8 



At) . At) _ .i (t) {2Ma) 



m 



m 



{E,v 2 (t)} = - 



_9_ _7_ 
16 + 16 



i 15 1 

cos(u;t) mu 2 v^(t) H sin(u;t) to v 2 (t) ^ — -mu 2 v l {t) = v 2 (t) 



(2.506) 

In spite of the fact that the physical energy E does not generate time evolution, it is 
conserved: 

dE _ dE . x dE . 2 dE 



dt -dv^ + ^ + -dt = ° 



(2.51) 
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This potential for a disagreement between the physical energy and the Hamiltonian which 
generates time evolution is a general feature of reduced canonical systems. 



3. Correspondence With Gauge Theories 

The purpose of this section is to study a simple model in which the various features of 
the construction and the context in which it takes place can be clearly understood. The 
model is scalar QED in flat space, the Lagrangian for which is: 

C=- 1 -F ilv + (fy - ze A^j 0* (> + ie (3.1) 

Here F^ v = d^A v — d v A^ is the field strength tensor, <fi is a complex scalar field, e is the 
electromagnetic coupling constant, and we are using a spacelike metric. The Lagrangian 
is invariant under a gauge transformation: 

A^x) ^ A^x) - dfjfiix) (3.2a) 

4>{x) \— > exp ie6(x) 4>(x) (3.26) 

which is parametrized by an arbitrary real scalar function, 9(x). By imposing temporal 
gauge: 

A (x) = (3.3) 

one obtains a constrained canonical system which still possesses invariance under gauge 
transformations parametrized by time independent functions 9(x). We will show how this 
system can be reduced to a completely gauge fixed system containing only dynamical 
degrees of freedom, and we will explain why the resulting Poisson brackets can fail to be 
canonical. We will then change variables to obtain canonical Poisson brackets and infer 
the Hamiltonian which generates time evolution for the reduced system. 

It is a simple matter to show that in temporal gauge the momenta canonically conjugate 
to A}, 4> and 4>* are, respectively: 

E i (t,x)=A i (t,x) (3.4a) 
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ic(t,x) = j>*(t,x) (3.46) 
w*(t,x) = j>(t,x) (3.4c) 
What it means to be "canonical" is that the non-zero equal time Poisson brackets are: 

{A t (t 1 x) 1 E J (t,y)}=5 lJ 5 3 (x-y) (3.5a) 

{(P(t,x),n(t,y)}=5 3 (x-y) (3.56) 

[<P*(t,x),n*(t,y)}=5 3 (x-y) (3.5c) 

The Hamiltonian is: 

H = Jd 3 x[^E i E i + ^ Fij Fij + tt* tt + (d, - ie A^j ^ (<9; + ie A^j </>} (3.6) 

It is straightforward to check that this functional generates time evolution. That is, the 
Poisson brackets of the Hamiltonian with the canonical coordinates give relations (3.4a) 
through (3.4c), and the Poisson brackets of the conjugate momenta give the canonical 
formulation of the dynamical Euler-Lagrange equations: 

E { = H | = dj Fji + ie (f>* (d{ + ie A^j <f> - ie cj> (d { - ie A^j <f>* (3.7a) 

7r={ir,H}=(di- ie A^j (d, - ie A^j 0* (3.76) 

n* = {tt, #} = (di + ie A^j [di + ie A^j </> (3.7c) 

The non-dynamical Euler-Lagrange equation — the one obtained by variation with respect 
to Aq — is not realized through the definition of time evolution. It must be imposed as a 
constraint: 

diEi + ie (7r0-7r*0*) = (3.8) 

Of course the left hand side of this constraint is also the generator of infinitesimal, time- 
independent gauge transformations. 
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Although this system possesses a time independent gauge symmetry, it is still complete 
in the sense that the equations of evolution uniquely determine the fields at any time in 
terms of those on an initial value surface. To economize the notation let us refer to a general 
field as tj) a (t,x) and an initial value configuration as i^ a (x). The evolution equations (3.4) 
and (3.7) uniquely determine the former in terms of the latter. Even though we cannot 
exhibit this relation as we could for the coupled oscillator system of the previous section 
— cf. relations (2.34) — we can still represent it functionally: 

V>a(t,£) = tf[?l (t,x) (3.9) 
L J a 

Just as with the coupled oscillator system discussed in the previous section, it is the ^ a (x)'s 
that represent the dynamical degrees of freedom of the unconstrained system. Only their 
Poisson brackets are independent, for example: 



6V 



U a (t,x),Mt,y)}= d 3 u d\ 1 J ; U c {u)Mv)\ (3.io) 

The equal time bracket algebra remains canonical because of the way time evolution acts; 
in our notation, because of the special way the functionals \l/ ip (t, x) depend upon the 

L J a 

initial value configurations. 

The constraint (3.8) represents a relation between the i^ a (x)'s. The theory is reduced 
by identifying a gauge condition on the ^> a (a?)'s that can be imposed by a unique, field 
dependent transformation of the residual symmetry group. Together the constraint and 
this gauge condition serve to eliminate a conjugate pair of the ^ a (x)'s. This does not 
change the way time evolution makes the fields depend upon their initial configurations, 
it only fixes the values of a conjugate pair of these initial configurations. The result is 
to change the equal time bracket algebra (3.10) since the Poisson brackets of the initial 
configurations are replaced with Dirac brackets. 

After reduction the full set of tp a (t, ^O's provides an overcomplete description of the 
system at any time. Although it is always possible to identify a minimal set of reduced 
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(t,x) ( ^ ^ , 5^> V (t,y) 



variables whose equal time bracket algebra is canonical, making such an identification is 
not necessarily easy and we often settle for a set of reduced variables whose bracket algebra 
is not degenerate but also not canonical. (If the bracket algebra becomes degenerate it 
means that the reduced variables do not provide a complete description of physics. That 
is, we cannot use the reduced set of ipaS to recover the reduced ipaS which are the true 
dynamical degrees of freedom.) This set is typically the evolution of whatever initial 
configurations remain after reduction. For example, in reducing the coupled oscillator 
system of the last section we eliminated q2 and p2, leaving q\ and p\. The natural first 
choice for reduced variables was accordingly q\{t) and p\{t). We then discovered that it 
would have been better to choose either q\{t) + q<i (t) and p\(t) + P2(t) or qi(t) — q2(t) and 
Pl(t) —P2(t) because these variables have a canonical equal time bracket algebra while the 
original choice does not. We shall now witness a similar phenomenon in reducing scalar 
electrodynamics. 

A natural gauge condition compatible with (3.8) is: 



From this and relation (3.4a) we infer the following evolution for the divergence of the 
vector potential: 



Of course (3.11) does not completely fix the gauge; it is still possible to make time indepen- 
dent, harmonic transformations. This freedom can be eliminated with a surface condition. 
A typical choice for the latter would be setting the normal component of A^(0, x) to zero 
on the "surface at infinity," which could realized as the limit of successively larger spheres. 
Rather than burdening the formalism with a cumbersome limiting procedure we will in- 
stead subsume the surface condition into the asymptotic fall-off usually assumed for the 



d l A l (0,x) = 



(3.11) 




(3.12) 
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gauge invariant part of the vector potential in order to make the magnetic field energy 
finite. 

There is a very close relation between fixing the harmonic gauge freedom and dividing 
the vector potential — and all other vector fields — into "transverse" and "longitudinal" 
components. With the just-stated convention the transverse and longitudinal parts of a 
vector field fi(t,x) which obeys the asymptotic fall-off condition are is defined as follows: 

fT(t, x) = fi(t, *)--!^J d 3 yG(x; y) djfj(t, y) (3.13a) 
flit, ^) = ^ij d 3 y G(x; y) djfj(t, y) (3.136) 



where the Green's function is, G(x;y) = —(4it\\x — y\ 




Since the longitudinal components of the electric field and the vector potential are 
entirely constrained it is natural to take the transverse components to be among the 
reduced dynamical variables, along with the scalar field and its conjugate momentum. 
Although the Dirac brackets of these variables are initially canonical: 

{Aj (0, x), Ej(0, y)) = Sij S\x-y) + ^-^- G(x; y) (3.14a) 

{^(0,x),7r(0,y)} = 5 3 (x-y) (3.146) 

{^(0,x),n*(0,y)} = 5 3 (x-y) (3.14c) 

time evolution makes the equal time bracket algebra non-canonical. We cannot exhibit the 
relation in closed form the way we did for the coupled oscillator system of the previous 
section but we can evaluate the second time derivative at t = easily enough. This is 
because reduction affects only the initial value fields; the time evolution equations are 
unchanged from (3.4) and (3.7). We can therefore use the evolution equations to reduce 
any number of time derivatives of an initial bracket to brackets of the initial value fields, 
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which can then be evaluated using (3.14). For example, we compute: 



= ^ x), n(t,y)} + [<j>{t, x), (di - ie A,) 2 cj>*(t, |^ 



(3.15a) 




0*(O, 



(3.156) 

y)} 




m[—-teA[(y)]^(y) 




(3.15c) 



Of course the failure of the 4>-tt bracket to remain canonical is due to the fact that reduction 
makes the longitudinal components of the vector potential and the electric field depend 
upon the charge fields. 

One consequence of (3.15) is that the bracket matrix contains explicit time dependence 
in addition to being non-canonical. To see this note that the bracket matrix must be time 
dependent because the second derivative of one of its components fails to vanish. If this 
time dependence were exclusively implicit — that is, if it derived only from the bracket 
matrix's dependence upon the time dependent dynamical variables - - then the initial 
matrix elements would depend upon the reduced dynamical variables. Because they do not 
we infer that the bracket matrix must harbor explicit as well as implicit time dependence. 

We saw in the discussion associated with (2.8) that explicit time dependence in the 
bracket matrix precludes the existence of a Hamiltonian which generates time evolution 
for the reduced dynamical variables. This is a little strange because (3.6) is the conserved, 
gauge invariant energy functional for scalar electrodynamics. However, it is easy to verify 
that the change reduction effects in the longitudinal electric field prevents even the initial 
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brackets from agreeing with (3.4) and (3.7): 



{0(O,X),#}=7T*(O,X) 



(3.16a) 



+ e 





} 



r d 




-dx l 



(3.166) 



Note again that our result is stronger than just that (3.6) fails to generate time evolu- 
tion: there does not exist a Hamiltonian which generates time evolution for these reduced 
variables. 

As was explained in the last section, the problem with finding a Hamiltonian that 
generates time evolution derives from an unfortunate choice of the reduced variables. We 
emphasize that there is nothing dynamically wrong with the choice we made. It is the 
natural one, and it does provide a complete and minimal description of the physics of 
scalar electrodynamics. Our choice is not even particularly inconvenient for computation; 
we will see in section 6 that it can be written in the functional formalism almost as simply 
as any other choice. The only problem comes if we insist on an explicit operator formalism 
in which time evolution is generated by a Hamiltonian. The cure for this problem is a 
necessarily time dependent field redefinition to canonical variables. 

There are many sets of canonical variables but the simplest is surely that which is 
related to our set of reduced variables through a time dependent gauge transformation. 
The gauge parameter is: 



We will adopt the convention that the Greek or Latin letter which denotes one of the 
old variables goes over in the new variables to the corresponding Latin or Greek letter 




(3.17a) 



ao(t, x) = ie 



I 



d 3 y G(x; y) { tt (t , y) cf>(t , y) - tt* (t , y) </>* (t, $ } 



(3.176) 
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respectively. The new variables are: 



aj(t,x) = Af(t,x) 
eJ(t,x) = Ef(t,x) 
f(t,x) = eKp^ie9(t,x) (j>(t,x) 
p(t, x) = exp — ie 9(t, a?) J |tt(£, x) — ie ao(t, x) (f>*(t, a?) j 
f*(t,x) = exp ^—ie 9(t,x) <fi* (t, x) 



(3.18a) 

(3.186) 
(3.18c) 
(3.18d) 
(3.18e) 



p*(t, x) = exp ie 9(t, x) |tt*(£, a?) + ze ao(t, x) (f>(t, x) j (3.18/) 

Since the quantity cto(£, is a gauge invariant it has the same form in terms of / and p 
as (j) and tv: 

a (t,x)=ie J d 3 yG(x;y) [p(t,y) f(t,y) - p*(t,y) f(t,y)} (3.19) 

By differentiating these relations and then using (3.4) and (3.7) we infer the following 
evolution equations: 



af(t,x) = ef (t,x) 



T i 



H (t, x) = 



dxi dxi 



ctj (t, x) + ie f*(t, x) di + iea^tjx) f(t,x) 



— * T — * ^ — * 

— ief(t,x) di — iea^tjX) f (t,x) 
/(£, x) = w*(t, x) - ie a (t, x) /(£, x) 



(3.20a) 
(3.206) 
(3.20c) 



p(t, x) = ie a (t, x) p(t, x) + -Q^- iea i (*» 



_d_ 

dx l 



-ieaf(t,x) f*(t,x) (3.20a 7 ) 



f*(t,x) = 7r(t,x) +iea (t,x) f*(t,x) 



(3.20e) 



r d 

p*(t, x) = —ie ao(t, x) p*(t, x) + — - + ie af(t, x) 

lax 1 



_d_ 

dx l 



T + ieaf{t,x) f(t,x) (3.20/) 



Of course these relations are generated by the Hamiltonian (3.6), which in the new variables 
has the form: 

H = J d 3 x j^ef ef + ^ d i a j d z aj + ^a 9ja + P* P+ (d z -ieaf^ /* (di + ieaf^j f} 

(3.21) 
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The reason the old Hamiltonian can generate evolution for these variables and not for our 
first choice of reduced dynamical variables is that the field redefinition contains explicit 
time dependence through the time integration in (3.17a). 

It should now be apparent — and it can be, with some difficulty, checked using relations 
(3.18) — that the bracket algebra is canonical. That is, the non-zero equal time brackets 
are: 

{af(t,x),ej(t,$)} =S lJ d 3 (x-y) + — — G(x;y) (3.22a) 

{f(t,x),p(t,y)}=5 3 (x-y) (3.226) 

{f^x)^{t,y))=5\x-y) (3.22c) 

It should also be apparent that the canonical formalism we have constructed is just that 
which follows from the invariant action (3.1) by the imposition of Coulomb gauge: 

diAi(t,x) = (3.23) 

Relation (3.19) comes from solving the constraint equation: 

didiA (t, x) = ie j>*(t, x) (f>(t, x) - <p(t, x) <f>*{t, x) (3.24) 

subject to a surface condition at spatial infinity which fixes the freedom to perform time 
dependent, harmonic gauge transformations. 



4. Perturbative Gravity Around Flat Space On T 6 x R 

This section is divided into four parts. In the first we describe the canonical formalism 
for gravity in a general closed spatial manifold. The second part introduces the mode 
and tensor decompositions we shall use for T 3 x R. In the third part we apply this mode 
decomposition to perturbation theory around flat space. It is here that we impose the 
constraints and fix the gauge to obtain the reduced theory. In the final part we obtain 
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a reduced canonical formalism and we show that our result agrees in the limit of infinite 
toroidal radius and localized initial value data with that obtained by A.D.M. [3] for open, 
asymptotically flat space. 

(1) — Description of Canonical Formalism. 

We define the lapse N° and the shift N l via the invariant interval: 

ds 2 = - (AT ) dt 2 + 1%3 [dx l + N l dt^ (dx j + N J dt^J . (4.1) 

This implies the 4-metric g^ v and its inverse g^ v are: 

- (N ) 2 + N k N l ~f k , N k ~f kl 
9nu=\ h 3 | (4-2) 



-1 



N3 
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[IV 



The usual Hilbert action for gravity 



(4.3) 



S = J d*x 
can be written in canonical form as: 

S = [ d 4 x 



(4.4a) 



(4.46) 



An integration by parts was used to arrive at (4.46) from (4.4a) and the following definitions 
were used: 



^ (JO*. 



= I V " v - V ' ) ('/,/ " -V;,./ - .V/.;,) (4.5a) 



K 



(4.56) 



Hi = -2 7lJ jv Jl . r (4.5c) 
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In the previous expressions a semicolon indicates covariant differentiation on the spatial 

sections using the connection compatible with the 3-metric, 7^; 1Z is the Ricci scalar 
formed from 7^. 

In these variables the Hamiltonian is, 

H = J (fixN^H^ (4.6) 
Variations of it with respect to n lJ and 7^ give us the evolution equations, 

iij = ^ N ° [*u - \-rij*) + N i;j + N r,i ( 4 - 7a ) 



fiU = - ( n ij _ lfj n ) + ^-N ^ (n lm n lm - J n 2 



K 



2 ^- 2 1 -j ■ 2 ^7 i y "im 2 



,-k 1 (4.76) 

while variation with respect to gives the constraint equation, 

T~Cfi = (4.7c) 

Of course (4.7a) is just a restatement of the definition (4.5a) of the conjugate momen- 
tum. Relations (4.76) are canonical versions of the six g^ Euler-Lagrange equations; the 
constraints (4.7c) are linear combinations of the four g^Q equations. 

We imagine the volume gauge to have been fixed by specifying the lapse and shift, 
possibly as functionals of the 3-metric and its conjugate momentum. Such a gauge condi- 
tion eliminates the ability to perform diffeomorphisms which are locally time dependent, 
as witness the fact that the Cauchy problem has a unique solution for fixed (and non- 
degenerate) lapse and shift. Just as with temporal gauge in scalar electrodynamics, our 
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gravitational gauge leaves a residual symmetry of transformations which are completely 
characterized by their action on the initial value surface, and by the condition that they 
do not affect the lapse and shift. 

Suppose we represent a general infinitesimal diffeomorphism, i— > + 9^{x), using 
the parameter 9^{x). It is a simple exercise to show that the 4-metric is changed by the 
following amount: 

Se (<T ) = O^p 9 pv + 9% - g^ p 9? . (4.8a) 
h = 9 P v ~ 0% g» P - gpv, P 9 p . (4.86) 

By requiring that Ogg^ = we see that the residual transformations, 6 9 (t,x), are 
characterized by their initial values, 9^(x), and by the following evolution equations: 

N° 

0° = 6°; N* -j^e p (4.9a) 



6* = 6* jNi-eVj^ 1 (n ^ +-JP-QPN 1 -N\ p QP . (4.96) 

As with temporal gauge scalar electromagnetism, the constraints generate residual sym- 
metry transformations. That is, if we define, 

ftp] = Jd 3 x [9°(x) A^(0,f)7^(0, x) + 9\x)H i (0 1 f)} , (4.10) 

then explicit calculation shows that 

{%, n[9]} = 9 \ % 3 + 9 k ;j % k + % 9° + 9% 7 jfc N k + 9 ^ % k N k (4.11a) 

=-V]N (4 - 116) 

{tt^', H [fl] } = n j 0° + ^(0 N° + 9 °<* N°>J + 9 ^ iV°' z ) 

(e >\ N° + 29 °' fc N° >k j - (iV l 5r fcj + N J 7? ifc ) fc 
+ 7? y iV fc k - 9 l . k 7T kj - 9 J . k 9 lk + (9 V 9 k ). k (4.12a) 

-- W*9- (4 - 126) 
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Note that this is not a definition. The left hand sides of (4.11) and (4.12) are defined 
by (4.10), (4.56) and (4.5c), while the identifications on the right hand side are made by 
applying (4.8) to the canonical coordinates and taking any time derivatives from (4.9). 

(2) — Mode Decomposition on T 3 x R. 

Now that the canonical formalism for a general space with closed spatial sections has 
been described it will be specialized to the treatment of T 3 x R. The coordinate ranges 
are t e R and < x % < L. The points x l = and x l = L are identified. Any function 
f(t,x) can be decomposed in modes in the following way: 

oo oo oo _ 

f(t,x) = L- 3 / 2 E E ^p\i^n ■ x] f(t,n) (4.13a) 



m=— oo ri2=— oo ri3=— oo 
-L r L r L 



f(t,n) = (2tt) -3 L~ 3 / 2 f dxi f dx 2 [ dx 3 exp f-i^n • xl f(t,x ) (4.136) 

JO JO JO L L J 

Note that when /(£, x) is real we have f*(t,n) = f(t, — ft). 

In representing tensors such as 7^ and 7r* ? it is convenient to decompose the index 
structure in a way that depends on the mode number. Let us define the 3-momentum, the 
transverse projection operator and the longitudinal inversion operator as follows: 

k = —n (4.14a) 

Ij 

Tij ee 8^ - (4.146) 

L ij = tij - 1^ ( 4 -!4c) 

For 7^ one can decompose any symmetric 2-tensor into three component pieces: 

fjk = f/k + fjk + i(fj h + Ik kj) (4.15a) 
If ee ^T lk T 3l - ^T fc ^ / w (4.156) 
/J ee l - Tij T kl J kl = I Tij J* (4.15c) 
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f.i = -p L .ikfjk (4.15d) 

— * 

Note that for each ^ there are two independent transverse traceless components 
j^i three longitudinal components /j, and one independent transverse component / *. 
Of course for k = all components satisfy the transversality condition. We therefore 
decompose the zero mode tensor into five transverse traceless components and one trace: 

1 



(4.16) 



We can carry the decomposition over into position space through the inverse transform as 
follows: 

fij = 3^ r ^*i + fij + fij + {fi,j + fj,i) 
f tr (t) = L-V 2 f tr (t) 



//j(t,f) = L- 3 / 2 X>p 
n 

4(t,f) = L- 3 / 2 ^exp 



i—n ■ x 



r 2tt_ 
%—n ■ x 

Lj 



■ 2tt ^ _ 
%—n ■ x 



7ij(t,n) 
fi(t,n) 



(4.17a) 
(4.176) 
(4.17c) 



(4.17d) 
(4.17e) 



/ 2 (t,x) = L- 3 / 2 ]Texp 

Note that the longitudinal and transverse components contain no spatial zero modes while 
f tr is all zero mode. The transverse traceless components alone contain both zero and 
non-zero modes. 

(3) — Perturbing Around Flat Space 

Since n^ v is a solution of Einstein's equations in T x R we can perturb around flat 
space, g^y = rj^ v + nh^ v . (We define the constant n 2 = 16nG.) The corresponding 
expansions for the various canonical variables are: 



lij 



— "ij + K h{j 



TV 



K 



V 



(4.18a) 
(4.186) 
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N° =1 + nn (4.18c) 
AT* =o + nn 1 (4.18d) 

We refer to h ij , p 1 ^ , n° and collectively as the weak fields. By convention the background 
metric is used to raise and lower indices on the weak fields. Since the background metric in 
this case is r\^ v it is irrelevant whether the spatial indices of weak fields are up or down, and 
raising a temporal index merely flips the sign. Note that the placement of k's in (4.18a) 
and (4.186) implies that h{j and p l] have the same bracket or commutation relations as 
7y and tv ij . 

If we expand the equations of time evolution, (4.7a) and (4.76), and then segregate 
according to tensor components, the following equations result: 

hfj = 2p« + 0(k) (4.19a) 

Pfj = \ V 2 h% + 0{k) (4.196) 

h^O + Oin) (4.20a) 

pt = _2V 2 n° + O(k) (4.206) 

V 2 hi + hj^i = -2^ + V 2 Hi + rijji + 0{k) (4.21a) 

V 2 Pt +p j7]t = + O(K) (4.216) 
■ h tr = _ p tr + D ( K j (422a) 

p tr = + O(k) (4.226) 

In these relations we have implicitly regarded the various weak fields as being of order 
one. This is not really correct because not all the fields are independent. Even in a theory 
without local symmetries we could use the equations of time evolution to express the weak 
fields at any time as functionals of the initial weak fields. It is traditional in this case to 
develop perturbative solutions as though the initial value configurations are of order one 
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in the coupling constant. The scheme is more complicated in a theory which possesses 
local symmetries because then one must impose a volume gauge condition in order to 
define a canonical formalism. Further, the canonical formalism so obtained possesses a 
set of constraints upon the initial value configurations and also, typically, a local but 
time independent residual symmetry. This residual symmetry is fixed by imposing gauge 
conditions on the initial weak field configurations. In our case we shall find it convenient 
to imagine that the surface gauge conditions are of order one, but we shall allow for the 
possibility of higher order terms in the volume gauge conditions. The constraints are solved 
perturbatively on the initial value surface to express the initial values of the constrained 
fields as power series expansions in functionals of the initial values of the unconstrained 
fields, regarding the latter as of order one. One then solves the perturbative equations of 
time evolution as for a theory without constraints but remembering that not all the initial 
configurations are of order one, and that the volume gauge conditions may also supply 
higher order terms. 

The four constraints can be expanded as follows in powers of the weak fields: 

Ti,Q = —(ji^i — hijij^j + (^—hh fc — —hhj k j ~ h,jhjk ~ ^ij^ij^k h-kih-ijj hijhkjjij ^ 
+ (pijPij ~ \p 2 ) + (~h 4 hj + ^h^hijj + hi l] k h lhk - -h ijyk h kj ^ + C(«)(4.23a) 

Hi = ~Pij,j - 2 ( h ijPjk) k + hj kj iP jk + O(k) (4.236) 

Substitution of the tensor decomposition (4.17) reveals that the Ho constraint determines 
the weak field h 1 : 



V V = kQq h l \ p u ; h\ p l ; h, p; h tr , p tr (4.24a) 



M ^tt. ut ^t. -l tAt Jut 



Qo = (~2^,fc + 2^i fc >i + hjhjk + hijhij k - h ki hijj - hijh k j^j ^ + (-PijPij + -p 

(\ h ,i h ,i - \ h ,i h ij,] - \ h ij,kkj,k + \ h ij,khkj,i) + °( K ) (4-246) 



+ 
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Similarly, the Hi constraint gives an equation for the weak field pf 



V 2 ^ + pjji = nQi [h tt , p u ; h\ p l - h, p; h tr , p tr j 



(4.25a) 



Qi = - (hijPjk) k + 2 h jk,iPjk + °( K ) (4.256) 

We can solve perturbatively for h l and pi because these weak field components contain no 
zero modes and the Laplacian is therefore a negative definite operator. However, we must 
first subtract off the zero mode parts of the sources Q^. For any function /(£, x) we define 
its non-zero mode part as: 

-L pL r L 



/ NZ (t,x) = /(t,f)-L- 3 / d Vl f dy 2 [ dy 3 f(t, 

Jo Jo Jo 



y) 



(4.26) 



To solve for and pi one simply inverts the Laplacian on the non-zero mode sectors of 
(4.24a) and (4.25a): 

h*=^ Q^[h tt J\h\p\h^h tr J r ] (4.27a) 

Pi = ^ Lij Qf [h u , p u - h\ p*; h, p; h tr , p tr ] (4.276) 

and then substitutes the resulting equations to re-express any h t} s or p^s which appear in 
the sources. For example the first iteration gives: 



K 



•»NZ 



Ull „EI. __r) 1N|Zj r> L - h — T ■ ■ NZ • h tr r> tr 

H ,p ,^2^0 ' V > a i ^2 l 3 ' 



K T ^NZTrtt tt. K n NZ t, K t NZ. h tr tr 
Pi = =2 Lij Qj h ,p ,^2^0 jPj^j =2^kl^l i h ,P 



(4.28a) 



(4.286) 



^2^13^3 [ ,F ' V 2 U V 2 " 
Of course there are still a*'s and p^s inside the new sources — though space presents us 
from displaying it explicitly — but whereas these fields might appear at order k on the 
right hand side of (4.27) they cannot appear before order k 2 on the right hand side of 
(4.28). Because each iteration moves them to a higher order in k we can obtain in this 
way an asymptotic series solution as a functional of h u , p u , p t , hi, h tr and p tr . 
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Although we have just seen that the constraints completely determine h l and Pi it 
is not quite true that constraining and pi completely enforces the constraints. There 
remain the zero modes. One can see by direct integration that although the zero mode 
constraints are free of terms linear in the weak fields they are not trivial at the next order, 
even when h? and pi are set to their constrained values: 



d 3 xH = J d 3 xQ (4.29a) 



-l L3 (P tr f + J ^{44■ + ^4>4>} + «Co[/l < ^p ^^ ;/l,p < ;/l ^ ^p <r ] (4.29c) 

J d 3 H l = J d 3 x Q l (4.30a) 
= j d 3 x \pj k hj k i + l -p l h*i + PjJ h kM - 2p Jjk h jM | + 0(k) (4.306) 
— > f d 3 xp% h% . + K d \h tt ,p tt - hj- h tr ,p tr ] (4.30c) 

(The functionals h tt ,p u ;h,p t ;h tr ,p tr are of cubic order and higher in the remaining 
weak fields.) A consequence is that there are solutions to the linearized field equations 
which can not be perturbatively corrected to give asymptotic solutions to the full field 
equations. This phenomenon is known as linearization instability, and it afflicts grav- 
itational perturbation theory whenever the background possesses Killing vectors* on a 
spatially closed manifold [17,18,19]. 

The linearization instability is sometimes regarded as a non-trivial obstacle to the 
development of perturbation theory. This is not correct. We need merely to restrict to those 
linearized solutions which satisfy the first non-trivial parts of the four zero mode constraints 



* The number of constraints which lack linear terms is equal to the number of Killing vectors. 
We have four because only the four translations give global Killing vectors for flat space 
on T 3 x R. Lorentz transformations — which also give Killing vectors for flat space on 
R 3 x R — do not respect the identification of T 3 x R. 
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and then develop systematic corrections as usual. Because our strategy is different for the 
global Hamiltonian constraint (4.29) than for the global momentum constraint (4.30) we 
shall discuss them separately. 

At quadratic order in the remaining independent weak fields the global Hamiltonian 
constraint is the difference of two manifestly positive quantities. This means we can solve 
it explicitly as follows: 



p tr = ±V6L-V 2 



d 3 x (pfjpfj + ±h? Jik h? Jik ) + kCq 



h tt J\hJ;h tr J r 



(4.31) 



The issue of choosing the sign in (4.31) commands considerably more attention than it 
deserves. The constraint equation does not fix it, and either choice is allowed classically 
- the positive sign corresponds to a contracting universe while the negative sign gives an 
expanding universe. If we are to avoid imposing extraneous conditions, and especially if 
quantum mechanics is to recover classical results in the correspondence limit then we must 
include both signs. This is achieved by using a multi-component wavefunction(al): 

(4-32) 



The action of the the wholly constrained operator p tr on ^ is defined by (4.31) with the 
plus or minus sign respectively. A purely contracting universe would be represented by 
= whereas a purely expanding universe would have \I/ = 0. We shall see in section 6 
that the straightforward application of the Faddeev-Popov technique for gauge fixing re- 
sults in the absolute value of an operator which causes the inner product to segregate into 
a manifestly positive contribution from each component. * 

There is also the issue of perturbatively iterating (4.31) to achieve an asymptotic series 
solution which is free of dependence upon p tr . We first define the zeroth order energy: 



Eo= I d\\p%p% + \h%^ k \ (4.33) 



* We wish to suggest that the same procedure be applied whenever discrete choices must be 
made in solving constraints. 
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and then expand the square root: 

= ±V6L- 3 / 2 I ^/Wq + kSq [h tt ,p tt ; h, p*; h tr J r \ } (4-346) 

Assuming that it is fair to regard the ratio Cq/Eq as of first plus higher orders in the weak 
fields, we then obtain an asymptotic series solution by iteration. For example, the first 
iteration gives: 

p tr = ±V6L- 3 / 2 | ^E~q + kSq p tt ; h, p 1 ; h tr , ±V§L-*/ 2 {^/E~q + kSq^ J (4.35) 

As with and pi, successive iterations push to ever higher orders any dependence of the 
right hand side upon p tr . 

Of course the iterative solution for p tr will result in nonsense, even perturbatively, if 
Eq can be made to vanish without Cq vanishing at least as rapidly. It turns out that this 
cannot happen for three reasons. First, we will shortly see that hj and p^ can be gauged 
to zero. Second, Eq is a sum of squares of all the remaining variables except for h tr and 
the zero modes of h u . Finally, the dependence Cq inherits from Hq implies that each of its 
terms must vanish at least quadratically with p**- and/or the non-zero modes of hfy To 
see this last point note from substituting the weak field expansions (4.18) into (4.56) that 
Hq consists of terms quadratic in p^j with any number of /ijj's and other terms which are 
free of pjj but contain at least one differentiated hij. Upon integration over T 3 each of 
the pure h{j terms must contain at least two non-zero modes of h^. Constraining and 
Pi to zero can result in terms which have any even power of the remaining components 
of Pij but it cannot result in odd powers of the momentum nor can it introduce pure 
terms which fail to possess at least two non-zero modes. Upon gauging p l to zero we see 
that every term in Cq must either contain a positive even power of pf- or p tr ', or it must 
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contain at least two non-zero modes of hfj. It follows that whenever Eq vanishes Cq must 



vanish at least as rapidly, so the ratio Cq/Eq can be legitimately regarded as of order one 
and higher in the weak fields. 

The three global momentum constraints cannot be imposed this way because we see 
from (4.30c) that their quadratic parts are not differences of squares. Our strategy is there- 
fore to leave them as constraints upon the classical initial value data or, in the quantum 
theory, upon the space of states. We can get away with this for three reasons. First, their 
imposition is not necessary in order to construct a reduced canonical formalism with a non- 
zero Hamiltonian. This was obviously not the case for the global Hamiltonian constraint. 
Second, the global momentum constraints remove no negative energy modes, unlike the 
global Hamiltonian constraint. Finally, the symmetry generated by the global momentum 
constraints consists of constant spatial translations on T 3 . Since these form a compact 
group there is no need to gauge fix them in the functional formalism. Additional support 
for the viability of not attempting to reduce the global momentum constraints can be 
found in the close parallel with the N = N constraint of closed bosonic string theory in 
lightcone gauge [20] . 

We turn now to the issue of gauge fixing. Since we wish in the end to compare our 
results with those of A.D.M. [3] we shall of course need to follow them in the choice of 
gauge. Their perspective was slightly different from ours: whereas we impose the volume 
gauge by choosing the lapse and shift, and then fix (most of) the residual symmetry with 
gauge conditions on the initial value surface, A.D.M. impose a volume gauge condition on 
the weak fields h{j and p lJ and then use the evolution equations for the frozen variables 
to determine the lapse and the shift. We can obtain the same result by merely choosing 
our lapse and shift, and our auxiliary surface conditions, so as to agree with A.D.M. The 
distinction between the two methods is important only to Faddeev-Popov gauge fixing in 
the quantum theory which was developed years after A.D.M. wrote. In our notation the 
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conditions favored by A.D.M. [3] are: 1 * 



hi(t, x) = 



p\t,x) = 



(4.36a) 



(4.366) 



We shall accordingly begin by showing that the residual symmetry allows the perturbative 
imposition hj = and p t = 0. We then argue that nP(t, x) and n l (0, x) can be chosen so 
as to perturbatively enforce the A.D.M. condition (4.36). 

To properly organize the notion of a perturbative transformation we insert a factor of 
k into the infinitesimal transformation parameter, 9^ = kt^. By substituting this and the 
perturbative expansions (4.18), and by iterating (4.11) we obtain the following expressions 
for the non-infinitesimal but perturbatively small transformations of h^j and Pij: 

\j — hij = — (ry + Tjj) 
1 

2 ( 



T ,i "./ j 
= ~( T i,j + T j,i) + ftWij 
^0 



- f J n % + ?° AJ f° - 2py f + pf°) + 0(k 2 ) (4.37a) 
^,p tt ;^,p^,p;^,p fr ;r°,f] (4.376) 



h t \p u -,h\p\hS-,h tr ,p 



tr ~tr~0 



■T ,T 



(4.38) 



We use the non-zero modes of % to perturbatively enforce = by iterating the equation: 

> - <r -° -1 (4.39) 



^ + ^2^j°k n jk 



~P\p tt ;h\p t ;hM tr ,p tr ^°,? 

The zero modes of are not fixed because they are conjugate to the global momentum 
constraints which are not being reduced. We use the non-zero modes of r° to perturbatively 
enforce p t = by iterating the equation: 

r° = -^- ^Tij Vf? [h U ,p U -h\ p *; ft, p; h tr , p tr ; f°, f] (4.40) 

* The component fields h{ and used by A.D.M. actually contain zero modes, unlike ours, 
and these zero modes were given non-zero values determined by the asymptotically flat 
boundary conditions [3]. 
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We will henceforth drop the prime and assume that hj = = p t . 

There remains the zero mode of f°. We shall use this to enforce h tr = although the 
argument for being able to impose this condition is more subtle. First, note that since the 
Wheeler-DeWitt symmetry must be gauge fixed [21], and since the subgroup of constant 
time translations is not compact, we do not have the option of declining to enforce some 
zero mode gauge condition. Second, note from (4.37a) that h tr is the quantity affected 
to lowest order by a constant time translation. Let us label such a transformation by the 
single parameter z: 

J-L-»/A?« W (4.4!) 

Whereas the parameter 9^ which enforced hj = = p l was of order k we need z to be 
of order one. Even so, the fact that p tr is constant to lowest order allows us to obtain a 
perturbative expression for the result of a non-infinitesimal shift: 

h tr ' -h tr = ^ r + K z\ i ti tt ,p tt ;h tr ,p tr -z\ (4.42) 

We see that the desired condition can be imposed formally by iterating the equation: 



z = -6 



ptr ^tr 



h tt ,p u ;h tr ,p tr ;z (4.43) 



We include the qualifier "formally" because the transformation is obviously singular when 
p tr vanishes. Of course p tr is not an independent degree of freedom; and we see from (4.31) 
that it is about as protected from vanishing as it is invar iantly possible to get. However, 
if all the modes of p^ and all the non-zero modes of hfj vanish then p tr vanishes, but we 
can have a non-zero h tr . In this case both hij(t, x) and Pij(t, x) are constant in space and 
time, and no temporal translation exists which will enforce h tr = 0. 

Our procedure is to go ahead and impose h tr = anyway. Considerable justification 
for this course derives from the close analogy to imposing the gauge g°(r) = t for a massless 
free particle whose position and momenta are Q M (r) and p^(r) respectively. In this system 
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the constrained variable, po(r) has an ambiguous sign which necessitates a 2-component 
wavefunction. Just as with gravity, the gauge condition conjugate to this constrained 
variable is singular for constant field configurations. This resolves itself in the quantum 
theory by the gauge fixed inner product acquiring a Faddeev-Popov determinant which 
endows the troublesome sector of configuration space with zero measure [21]:* 



1> 



5[q — t ) abs ( pq 



d 6 q { if; + *(t, x) <P + (t, x) - x) <P + (t, x) 



+ - / d*q< V~ (t,x)<f>~(t,x) *(t,x)<l>~(t,x) } (4.44) 



We will see at the end of section 6 that Faddeev-Popov gauge fixing endows gravity with 
the same sort of inner product. The result obtained for the free particle has such uni- 
versal acceptance that we shall henceforth ignore the completely analogous problem with 
imposing h tr = on constant field configurations. 

It remains to show that we can choose the lapse and shift so as to enforce the A.D.M. 
gauge conditions (4.36) for all time. To see this it suffices to apply the constraints to 
evolution equations (4.206), (4.21a) and (4.22a): 



p* = -2V 2 n° + k A/" [h tt , p tt ; n°, n; p\ h, h tr j 



(4.45a) 



7 * 9 r 

V hi + hjji = V rii + rij.ji + nNi 



h tt ,p tt ;n°,n;p\h,h tr j (4.456) 

k tr = T L~ 3 / 2 v^Eb + « T [h tt , p tt ; n°, n; p\ h, h tr j (4.45c) 

Since we already have p*(0,af) = = h{(0,x) we will have the A.D.M. condition (4.36) if 
the weak field lapse and shift are obtained by iterating the equations: 



_ K A r0 



n = 



/i tt ,/;n°,n;0,0,a <r 



K 



•NZ 



/i tt ,p tt ;n°,n;0,0,/i* r ] 



(4.46a) 



(4.466) 



* We have taken the liberty to correct an error that appeared in formula (23c) of [21]. 
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Note that these equations only determine the non-zero modes of the lapse and shift. We 
propose that the zero modes be left one and zero, respectively, to all orders. Of course while 
(4.45c) — and the initial condition, h tr (0,x) =0 — determines h tr (t,x), this component 
field does not vanish after t = 0. 

(4) — The Reduced Canonical Formalism and Its and Correspondence Limit. 

In the previous section we succeeded in reducing the theory to the point where only 
the transverse-traceless fields survive. It is convenient, however, to view the system that 
results when p tr and 9 are not yet reduced. Because special care must be given to the 
zero modes we will do this in /c-space: 



tt 



p.J = h U + K Ul] {h T \p u ) + h lr h* - \p tr p$ + 0( K ' Z ) (4.47a) 



,tt Jt 



1 



2 utrZ tt 



1 



±r~ tt 



12 



tt 



2W ij (h t \p u ) + lh tr p i f + ±p tr h i f 



+ 0(k 2 ) 



(4.476) 



The explicit forms of Uij and W^ ? - are: 



U, 



'J 



1 



(2tt) 3 L 3 / 2 V Jr 2 3 



nr 



d 3 x <! e 



ik-3 



-h" ! -' 



tt 



tt 



-h tt 



tt 



- h tt h U 



tt 



tt 



2 lm,n lr,m 2 ln,m lm,r 2 n K m rl,m' r 2 nl,m rm,l 

-h U h tt - -h U X7 2 h U - -h tt X7 2 h tt - 2n tt n tt ] X 
2 n lm n nr,lm 2 nl lr 2 nl lr nl lr J 



(4.48a) 



(27r) 3 L 3/2 V T<nT J' r 2 TljTnr 



d 6 x 



-ik-x I tt 



h u _l_ n u h u 

lu mr > frm ""mn 



nm '"mr 



(4.486) 



and the evolution equations for the zero modes are obtained by setting k = after changing 
Tij — > Sij in the expressions for Uij and Wjj. 

This system is canonical for the same reason that the A.D.M. system is: the surface and 
volume gauges have been chosen so that the variables conjugate to each of the constrained 
variables - - that is, p t for h l and hi for pi - - remain zero for all time. The act of 
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reducing p tr and h tr spoils canonicity because h tr does not remain zero after the initial 
time. The mechanism is the same as we found in section 3 where the evolution of a 
non-zero longitudinal vector potential in temporal gauge broke the canonicity of scalar 
electrodynamics. Because it is only the trace components which break canonicity we know 
that it will be restored if we can transform to variables (X, P) with the same evolution 
equations except for lack of dependence on p tr and h tr . This transformation will necessarily 
be time dependent and non-local. The time dependence arises because the transformation 
must give (h u ^p u ^ at t = so that at this time (X,P) obey the same commutation 
relations obeyed by (h u ,p u ); but it must deviate from this at later times since we wish 
to eliminate (h tr ,p tr ) from equations (4.47). The non-locality enters merely because the 
transformations must depend on (h tr ,p tr ) and these are non-local, as witness equation 
(4.31). 

It is trivial to check that to the order we are working the following transformations 
possess the properties mentioned above: 

For u> = \k\ 7^ 0: 

sin(a;t) 



Pi f cos(ut) + sm(ut) 



Y — h 1 



tt 



K< 





sin(cjt) 











while for uj = 0: 



h { f cos(cut) - —Pi™ sm(ut) 



^h tr hA (4.496) 



PiM = p#(o) + H - Y2 htr Pif (°) 



(4.50a) 



XijiO) = htfip) - K^tp tr [hifOO) - tp^)] - l^fc.«(0)} (4.506) 

Applying the transformations (4.49) and (4.50) to the evolution equations (4.47) results, 

as required, in equations independent of both h tr and p tr : 

1 

2 
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Pij = -^k 2 Xij + KlAij{X, P) + 0(k 2 ) 



(4.51a) 



X^ =2P ij + 2 K WiAX, P) + 0{k 2 ) 



(4.516) 



and again, the equations for the zero modes are obtained by setting k = in the above 
after changing Tjj — > Sij in equations (4.48). 

As previously mentioned, the variables X and P are canonical at t = since at this 
time they are simply equal to h u and p**. Furthermore, they will remain canonical at later 
times since the evolution equations (4.51) are just those of the A.D.M. variables. 

The Hamiltonian that generates equations (4.51) in terms of the x-space variables is: 



H = I d°x{ P i:j P i:j + - X ijtl X ijtl + k 



(4.52) 



there are two reasons for writing the Hamiltonian in terms of {X, P) as opposed of (X, P). 
The first one is merely that it is in terms of these variables that the form of H is the 
simplest. The second reason is that in this form it is obvious that the Hamiltonian density 
implied by (4.52) is local. However, in order to derive the evolution equations (4.51) 
while still treating the zero modes properly, one must work in momentum space. Here the 
non-zero bracket relations are: 



For k,q^0, 



X tJ (k),P lm (-q)\ = -8^ 



1 



TnT, 



1, 



il 1 jm r > 1 ij- L lm 



TaTim ) + I Tj m Tji T-ijTi m 



(4.53a) 



For k = q = 



X ij (0),P lm (0)\ = - 



(4.536) 



So far we have succeeded in reducing the theory and extracting the proper Hamiltonian. 
We will next prove that in the appropriate limit the Hamiltonian of equation (4.52) goes to 
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that obtained by following the A.D.M. procedure. The appropriate limit is that in which a 
configuration in T 3 x R goes to the same configuration in an open space with flat boundary 
conditions. Explicitly, the limit in which the two treatments agree is that in which we take 
L — > oo with localized initial perturbations. We refer to this limit as the open space limit. 
It should be obvious that the proof reduces to showing that in the open space limit both 
p tr and h tr vanish; since if this is the case (X,P) become just {h u ,p u ) respectively and 
the Hamiltonian (4.52) already has the correct form. 

Let us begin then by examining Eq as defined by equation (4.33): 



E = I d*xW? j p? j +--KX k hZ 



1 



tt u tt 



(4.33) 



ij p ij ~ r 4 ij,k b ij,k 

Note that Eq remains finite in the open space limit, even though the range of integration 
increases from [0, L) to (— oo, oo). The reason for this is that localized initial perturbations 
guarantees that the integrand above has finite support. 

Now let us inspect the evolution equation for h tr (equation (4.22)) together with the 
constraint equation for p tr (equation (4.31)): 



Itr 



h 



-p tr + KL 



-3 



p tr2 = 6L~ 3 { Eq + k I d 3 xA ^J 1 - h tr J r 



(4.54a) 



(4.546) 



Equations (4.54) are iterative relations for h tr and p tr in terms of the tt fields after the 
gauge has been fixed and the constraints for and pi have been enforced. Equation 
(4.54a) can be integrated (again iteratively) to give: 



h 



tr 



dr 



-p tr + nL- 3 



d 6 x{ A\ 



h tt ,p tt ;h tr ,p tr 



(4.55) 



to see that both h tr and p tr vanish in the open space limit we must examine L~ 3 J d?xA\ 
and L~ 3 J d?xAQ closer. Let us then explore the L dependance of each of these two terms 
separately: 
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Ai is of second order and higher in the fields; therefore the highest power of L in 
J cftxAi occurs when the integral acts on constants (since both h u and p u have finite 
support). In the open space limit we can then replace J cfixA\ with L?M.\ where M.\ 
is at least of second order in h tr and/or p tr . 

The form of Aq is at least of third order in the fields (remembering that Eq is inde- 
pendent of L in the limit). Similar considerations as those mentioned for the case of 
A\ reveal that in the open space limit J cftxAo can be replaced by L^A4q where A4q 
is of degree three and higher in h tr and/or p tr . 

We can then, in the open space limit, write equations (4.55) and (4.546) as: 

t 



= J dr j -p tr + kMi [h tr ,p tr j | (4.56a) 



h tr 



p tr2 = L- 3 E' + k M'q [h tr ,p tr ] (4.566) 

where the primes have been put in to absorb an irrelevant factor of 6 that would otherwise 
appear multiplying the right hand side of equation (4.566). At this point it should be 
obvious to the reader that the perturbative solutions to equations (4.56) are h tr = and 
p tr = 0. For those that still have some doubts let us take the L — > oo limit and re-write 
equations (4.56) as: 

t 

h tr = Jdr{-p tr + Ka nm (h tr Y \p tr ) m } (4.57a) 


p tr2 = Ka' n , ml (h tr Y (p tr ) m ' (4.576) 

with n + m > 2 and n' + m' > 3. Each successive iteration of equations (4.57) brings with 
it positive powers of k. Therefore, to any order in perturbation in powers of k both h tr 
and p tr vanish in the open space limit. Thereby proving the correspondence between our 
method and that of A.D.M. 
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We conclude this section by re-stating the result: Our method of reduction gives a 
precise meaning to time and, perhaps more importantly, this time evolution coincides in 
the appropriate limit with that obtained by A.D.M. for open space. 



5. Exact Results For Minisuper space 

In this section we examine two simple models. The first one is that of pure gravity 
with a cosmo logical constant. The second example deals with the case a gravity coupled 
to a scalar field. In both of these models we will truncate the form of the metric to allow 
only certain zero modes to be present so as to render the models exactly solvable. 

— Zero modes of gravity with a cosmological constant. 

For the first example we choose the form of the 3-metric to be such as to have two 
homogeneous degrees of freedom: 

ds 2 = - N 2 (r)dr 2 + 6 2 / 3 (r) [e a ^ dx 2 + e~ a ^ dy 2 + dz 2 ~\ (5.1) 



the action in canonical form is: 

S : 

with 1Z defined as: 



J dr p b b + p a d-NK. (5.2) 



<i? k2 2 2 A 3k 2 



The equations of motion before the reduction is implemented are obtained by varying 
the action of (5.2) with respect to p a , a, and b respectively: 

K 2 

a = N-g I p a (5.4a) 

Pa = (5.46) 
3k 2 

b = -N — bp b (5.4c) 
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while varying N results in the constraint 71 = 0: 



3k' 



k 2 2 2 A 



bp( = 



(5.4e) 



Since the constraint (5.4e) will be enforced by singling out p^ we wish to choose the 
volume gauge by simplifying the equation for 6 as much as possible. The obvious choice is: 



N = 



b\ Pb \ 



(5.5) 



Now that the volume gauge has been fixed we fix the constraint and use equation (5.4c) 
to fix b in the following manner*: 



Pb = 



16 A 4 9 

+ —pi 



3 k 4 



3 6 3 



(5.6a) 



3k z 
1 + — 



(5.66) 



The reduction is now complete. The equations of motion for the physical fields become: 

~i -1 



a = 



K 2 p a 

6(*) 3 



'16A 4 

+ - 2 



3 k 4 



3b(t) 



3 Pa 



(5.7a) 



Pa = 



(5.76) 



the above equations are integrable and the Hamiltonian can be obtained from them: 

3 k 2 



H 



16 A 

J^ 4 



+ 



3b(t) 



3 Pa 



(5.8) 



Before going to the next example we wish to clarify one point: The choice of N ^ I 
signifies that the time evolution implied by equation (5.8) is not that corresponding to 



* For simplicity of exposition we have made a definite choice for the sign of p^, that cor- 
responding to an expanding space. As in section 4 the wave function really consists of 
two components, one for each of the two signs. Note also that we fix the surface gauge 
condition by choosing 6(0) = 1. 
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time evolution in flat space (we will see this point more clearly in the next example when 
we take the limit k — > 0). N was chosen so as to make equation (5.4c) exactly solvable. It 
was by no means a unique choice; for example had we chosen: 

1 



N = 



\Pb\ 



(5.9) 



equation (5.4c) would still be easy to solve but equations (5.4a — b) would have a different 
form: 

I 2 e Pa / M M 
a = k (5.10) 



Pb 



Pa = 



where prime denotes differentiation with respect to the new time parameter t. The Hamil- 
tonian then would be: 

.2 



H t = 3 K * c -3/4« 2 t ' 16A 



+ 



,3 Pa 



(5.11) 



i \/ 3k 4 ' sb(ty 

This gauge dependance of the Hamiltonian should come as no surprise since changing how 
we gauge fix N changes what we mean by time, thereby changing what we mean by time 
evolution. In our next example we will show how despite this freedom we can make contact 
with the results obtained in a theory for which gravity is not dynamical (i.e., in the limit 
k — > 0) and N = 1 always. 

— Zero mode of gravity coupled to a scalar field 

We start by looking at the zero modes of a massive scalar field coupled to gravity and 
we allow the 3-metric to have only one degree of freedom: 



ds 2 = -N 2 (t) dr 2 + exp 



-a(r) 
3 V 1 



dx 



*2 



(5.12) 



The action in canonical form is: 



S= dr 



Tv<p + pd- NK 



(5.13) 
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where (tt , p) are the variables conjugate to (4> , a) respectively and: 



n = 



1 



1 2 2 

a p 



e~ a + 



1 



m 



2 j2 



(5.14) 



with a 2 = YI-kG. 



By varying this action we obtain the unreduced equations of motion: 

a = -a 2 Ne~ a p 
p = N 

= Ne~ a n 



1 2 1 2 2 
- 7r — - a p 

2 2 1 



e- a - N 



1 

- m 
2 



2 j.2 



and 



ir = -Nm 2 (f)e a 



while variation with respect to iV gives the constraint equation 



1 



-n 2 - - 
2 2 



1 2 2 

a p 



e~ a + 



1 2 1 2 
-mo 

2 ^ 







(5.15a) 
(5.156) 
(5.15c) 

(5.15rf) 
(5.15e) 



We must now select a volume gauge condition to impose. Setting N = 1 will not do 
because one would be left with the task of solving equation (5.15a) explicitly for a. A 
gauge choice that simplifies this task is: 



N = — 
\P\ 



(5.16) 



It is obvious that the above choice makes the job of solving equation (5.15a) a trivial one. 
Note; however, that iV does not approach 1 in the limit k — > 0. We will have to account 
for this when comparing our results to those obtained for flat space. 

Having fixed the volume gauge we now proceed to reduce the theory by enforcing the 
constraint (5.15e) and fixing the value of a at r = *: 

1 



p = \l n' z + m 2 4> 2 e 2a 

a 



a = a t 



(5.17a) 
(5.176) 



* We again chose the sign of the constrained variable to give increasing a for increasing r 
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where we have chosen a(0) = 0. 

The equations of motion for the remaining variables are, after implementing reduction: 

(5.18a) 



a//t 2 + m 2 4> 2 e 2a 

-a m 2 (pe 2a , rio ,, 

7t = y (5.186) 
\Ar 2 + m 2 (/> 2 e 2a 

The gauge choice (5. 16) leaves the variables canonical. The Hamiltonian obtained by 
integrating equations (5.18) is: 



H T = a\Jit 2 + m 2 (/) 2 e 2a2T (5.19) 

the subscript r is there to remind us that this Hamiltonian describes evolution with respect 
to r which in the limit a — > does not go to the flat space t simply because iV does not 
go to 1 in that limit. 

To see how to recover the flat space result in the limit a — > 0. Let us examine the r 
evolution of <fi: 

d<j) 8H T 



dr dn 



(5.20a) 



a e 



a 2 Ht) 



> n 2 + m 2 02 e 2a2r(t) 
which in the limit a — > can be written as: 



d A = ^- (a a/tt 2 + m 2 2 e 2 «Mt)^ (5.20c) 
ot on \ * J 



., 7T 2 + m 2 (/> 2 — \/ 7r 2 + m 2 </> 2 (5.21a) 

Ot * OTT v 



9 2 a. 1 2 ^ 
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(5.216) 



which of course is the correct limit. 

We end this section by pointing out what we hoped to accomplish with these two exam- 
ples. The first example was meant as a simple illustration of the method for an admittedly 
simple model. In it we made evident the fact that the form of the reduced Hamiltonian 
rests on the choice of the lapse function (i.e., the Hamiltonian is gauge dependent). The 
second example was used to show how the paradox of dynamics is resolved by reduction. 
It was used to show that the limit of k — > gives the same results obtained for flat space. 

6. Correspondence With The Functional Formalism 

The purpose of this section is to show that expectation values and matrix elements in 
the reduced canonical theory can be expressed very simply in terms of the naive functional 
formalism of the unconstrained theory. The key to obtaining this result is that reduction 
affects only the allowed initial values of Heisenberg operators, not their subsequent time 
evolutions. We can therefore enforce reduction by gauge fixing on the initial value surface 
and use the unconstrained Hamiltonian to implement time evolution. Since evolution 
is carried out with the known Hamiltonian of the unconstrained theory we get the usual 
functional formalism and we need never go to the trouble of actually constructing a reduced 
Hamiltonian or inferring how the original dynamical variables depend upon the reduced 
canonical variables. We first derive the result for a general constrained canonical formalism, 
then we explain how it applies to the coupled oscillator model of section 2, to scalar QED 
in temporal gauge, and to quantum general relativity with fixed lapse and shift. 

We require first a notation for the unconstrained canonical formalism whose reduction 
leads to the z/(t)'s of section 2. In a gauge theory this would be the canonical formalism 
which results from the imposition of a volume gauge condition, that is, one set of conditions 
for each spacetime point. Fixing the volume gauge is enough to make time evolution unique 
but it typically leaves a residual symmetry of gauge transformations of the initial value 
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data. These residual symmetry transformations are generated by constraints which restrict 
the allowed initial value data. It is the act of fixing the residual symmetry and imposing 
the constraints which leads to the reduced dynamical system characterized by the ?/(£)'s. 

Let us consider the unconstrained system to be described by coordinates x a (i) and 
momenta 7r a (t) where a = 1, 2, . . . , L, and L = N + K . The variables are canonical, that 
is, the only non-zero equal time commutators are: 



x a (t),^(t) 



iS a , 



(6.1) 



Time evolution is generated by a possibly time dependent Hamiltonian, H(x, 7r, t): 



± a {t) = -i x a (t),H(x(t),n(t),t^) (6.2a) 

n a (t) = -i [ir a (t),H(x(t), tt(£), *)] (6.26) 

Equations of this form serve to determine the variables at any time in terms of the initial 
values, x a and n a : 

x a {t)=X a (x^,tj (6.3a) 

7T a (t) = II a (x, 7?, t^j (6.36) 

Finally, we assume the system is subject to a set of K constraints: 



C k (x{t)^{t),t)=Q 



k = 1,2, ...,K 



(6.4) 



To accommodate the coupled oscillator model of section 2 we have broken the usual con- 
vention of forbidding explicit time dependence in the constraints. We do require, however, 
that if Cfc (x, 7r, j =0 for all k then Cy. (x(t), 7r(£), t j =0 as well. 

States in the unconstrained formalism can be represented by their wavefunctions in 
the basis of position eigenkets at some fixed time: 



iM) = f d L i v(0 



(6.5a) 
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x a (t) 



(6.56) 



Note that these are Heisenberg states so they do not evolve in time, despite the fact that 
the wavefunction might be given in terms of the eigenkets of the position operator at any 
time. The inner product between two such states is the usual one: 



(6.6) 



The Heisenberg evolution operator: 

U(t 2 ,h) = rjexp ijT * dtH(x(t),ir(t),t) j (6.7a) 

(the symbol T denotes the ordering convention in which canonical operators at later times 
appear to the left of those at earlier times; coordinates stand to the right of momenta at 
equal times) enters when we wish to study operators at different times: 

x a (t 2 ) = U(t 2 ,t 1 )x a (t 1 )u\t 2 ,t 1 ) (6.76) 

or when we wish to specify the wavefunctions in terms of the eigenkets of operators at 
different times: 



(6.7c) 



t;t 2 ) = U(t 2 ,h) 

In these cases it is convenient to subsume the evolution operator into a functional inte- 
gral. Suppose we wish to study some functional 0[x, ir] of the canonical operators defined 
between times t\ and t 2 > t\. We can take the matrix element of its time-order product 
between states specified at t 2 and t\ using the following formula:* 



^25*2 



t(o[ 



X, TV 



x exp 



= J[dx(*) dTv(t') il>Z(x(t 2 j)0 

*2> t > t l t 2 >t'>t 1 

dt[n a (t)x a (t) - H(x(t),ir(t),t)} Vl(^l)) (6-8) 



* Functional integrals are sufficiently familiar in particle physics that we shall omit spec- 
ification of the detailed skeletonization which would be necessary to make our formulae 
well-defined [22]. This could only be done for quantum mechanics and free field theory 
in any case. Complete rigor has not yet been obtained for interacting quantum field the- 
ories in four dimensions, although it is completely obvious how to proceed in regulated 
perturbation theory. 
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Note that the symbols u x a (t) n and u n a (t) n stand for operators on the left hand side of 
this equation while they are C-number functions in the functional integral to the right. If 
the Hamiltonian is quadratic in the momenta then we can perform the Gaussian integra- 
tions over 7r a (t) and pass from the phase space functional formalism to the more familiar 
configuration space version. 

If we wish to specify both states at the same time — say t\ — then it is necessary to 
first evolve forward to some arbitrary point beyond the final observation in 0[x, 7r], and 
then evolve backwards to t\. The necessary formalism was worked out three decades ago 
by Schwinger [23] and has been studied more recently by Jordan [24].** If we denote the 
fields which implement forward evolution with a superscript or subscript "+," and those 
which generate backward evolution with a superscript or subscript "— ," then the relevant 
functional integral is: 



V>2;£l 



t(0[x,tt] 



X 



V> 2 *(x-(ti)) 



exp 



—i 



dt 





x+,n + 






x O 


exp 


V 











dx-(t)^dn (t') dx+(t) dn + (t') s(x-(t 2 ) - x+(t 2 ) 

*2>*>*1 t 2 >t'>t 1 t 2 >t>t 1 j 2>t '> tl 

{ n a (t) x a _ (t) - H (x- (t) , 7T - (t) , t) } (6.9) 

dt {n+(t) x%(t) - H(x+(t),7r+(t), t) } Vl (^(*l)) 

Note that we must use "+" variables in the functional O if the functional integral is to 
give the expectation value of the time-ordered product; using the "— " variables would give 
the anti-time-ordered product. 

The unconstrained matrix elements and expectation values we have described are in- 
adequate in two ways: they are typically divergent for the only interesting states, and they 
include information we don't want to know about unphysical operators. The first problem 
arises because one enforces the constraints (6.4) by requiring the states to be annihilated by 



** Although Schwinger and Jordan assumed that the initial and final states were free vacuum 
and that t\ was in the asymptotic past, the generalization to arbitrary states and times is 
trivial. 
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the Cfc's. The inner product (6.6) between two such states can be divergent on account of 
the integration over residual gauge transformations. This is even true when, as in the case 
for gravity, the residual gauge transformation of a coordinate x a involves the momentum 
up [21]. The second problem arises because it is really operators in the reduced canonical 
formalism whose matrix elements and expectation values we wish to study. These reduced 
operators have the same time evolution as the unconstrained ones, but they depend upon 
2K fewer initial value operators. 

It turns out that by implementing reduction we can also solve the problem with the 
inner product. The unconstrained theory is reduced by identifying K residual gauge con- 
ditions on an initial value surface which we take to be t = for definiteness: 

G k (x,7?)=Q (6.10) 

These conditions are arbitrary except for the requirement that the Faddeev-Popov matrix 
formed from commutation with the constraints: 



' (x, tt) = -i C k (x, 7T, o) , G e (x, tt) (6.11) 

should have a non-zero determinant. 

We can use the conditions G k = and Ci = to decompose the 2L operators of the 
unconstrained formalism into two commuting sets of canonical variables: 

{x,^} — [(q,p); (c,g)} (6.12) 

The iV q a, s and iV p^s commute canonically and are of course the quantum versions of the 
reduced canonical variables whose construction was described in section 2. The K g fc, s 
form similar conjugate pairs with the K c^'s. The two sets of variables commute with one 
another so that the only non-zero commutators are: 



q a 3b 



= iS a b (6.13a) 
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(6.136) 



The g fc, s are pure gauge variables which have the property of vanishing with the gauge 
conditions: 



9 k 6 K 



G 



= 



(6.14) 



The C£S are constrained variables; the constraint equations determine them as functions 
K e(QiPj °f the reduced canonical variables. The Coulomb potential of electrodynamics 
and the Newtonian potential of gravitation are examples of constrained variables. When 
acting on invariant states — that is, upon states which are annihilated by the constraint 
operators — there is no difference between the c^'s and these functions: 



V'inv;*) = 



V'inv! t 



V'inv > t 



(6.15a) 
(6.156) 



Except for the fact that we have distinguished between the constrained variables and 
the constraints themselves — that is, we have allowed for the possibility that K£ ^ 0, 
and that the reduced canonical variables might not commute with the constraints — the 
decomposition (6.12) is a standard one in the theory of constrained quantization [25]. The 
reduced operators whose matrix elements and expectation values we really wish to study 
are: 

x?(t) = X a (x,n,t) ^ ^ (6.16a) 

V / g=0 , c=K 

7r r a {t) = U a (x,7r,t) ^ ^ (6.166) 

V / g=V , c=k 

Note that they have the same time evolution as the unconstrained operators (6.3). 

Now consider the inner product of two states which are annihilated by the constraints. 
Since the potential for divergences in the naive inner product comes from integrating over 
residual gauge transformations we can avoid the problem by surface gauge fixing. In the 
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case where the gauge conditions and the Faddeev-Popov determinant depend only upon 
the coordinates this is accomplished by formally inserting unity in the form: 



d K 9 exp 



G\X,7T 



abs< det 



M^i ( X, 7T 



exp 



(6.17) 



One then changes variables and drops the typically divergent integration over gauge param- 
eters. Although the question of operator ordering arises forcefully when either the gauge 
conditions or the Faddeev-Popov determinant are allowed to depend upon non-commuting 
operators, it is natural to expect that the same procedure works generally. This expecta- 
tion is fulfilled in several interesting cases, including that of perturbative quantum gravity 
[21]. We therefore propose the following reduced inner product: 



^25*2 



G 



abs | det 



M U X, 7T 



(6.18) 



where it is understood that the issue, if there is one, of operator ordering is to be resolved 
on a case-by-case basis. Note that in view of (6.17) this inner product is independent 
of the choice of gauge for states which are annihilated by the constraints. In fact it 
amounts to a realization in the unconstrained formalism of the standard inner product of 
the reduced theory [25]. Enforcing this correspondence is an important determinant of 
the ordering convention for the gauge fixing paraphernalia; another important condition is 
the persistence of Hermiticity when we drop the requirement that the wave functions are 
annihilated by the constraints. 

The relation we seek between invariant operators in the unconstrained theory and the 
corresponding reduced operators holds in the gauge fixed inner product (6.18). Although 
invariant operators generally depend upon the g fc, s and the c/s in addition to the reduced 
canonical variables, when such an invariant operator acts upon an invariant state we can 



use the constraints eg = Kg 



O 



mv 



q,p 



X, 7T 



V'inv ! ^ ) ^inv 
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X, 7T 



C=K 



(6.19a) 



o 



mv 



(6.196) 



Note that this is true even when the c/s are not themselves invariant. Once the c/s have 
been eliminated there is no obstacle to commuting the g fc, s through the q a, s and fys to 
act on the gauge fixing delta function. We therefore obtain the result that the expectation 
value or matrix element of an invariant operator in the presence of invariant states is the 
same as the expectation value or matrix element of the corresponding reduced operator: 



o 



mv 



X Tl 7T 



G 



(6.20) 



To reach the functional formalism we consider the action of the gauge fixing parapher- 
nalia on the ket to be a state in the unconstrained theory: 



G (x, 7T j abs < det M^g , 



(6.21) 



and apply the functional formulae (6.8) and (6.9). It is clearly convenien for this purpose 
to either assume that t' = or else to fix the residual gauge on the surface t = t' . The 
relevant formulae for matrix elements and expectation values are, respectively: 



^25*2 



T[0[x,n] 



x exp 



dx(t) dixit') i>l(x(t 2 ))0 

*2>*>*1 t 1 >t'>0 

i dt {yr a (t) x a {t) - H(x(t),ir(t),fj } 



x 6 



K 



(6.22a) 

7r(ti))] absj det [M fe ^x(ti), 7r(ti))] | 



^25*1 



t(o[. 



X, TV 



X 



x C 



x 5 



K 



oV - ^') ete + (£) d7r + (t') <y(x_(t 2 ) 

*2>*>*1 t 2 >t'>t 1 t 2 >t>t 1 t 2 >t'>t 1 

2 dt |tt- (t) x a (t) - H[x-{t), Tc-(t),t) } 
r {yr+(t) «4(f) - ff(z + (t), 7T+(t), f) } 
G(x(*i),7r(*i))] abs|det[M w (x(ti),7r(ti))] Ui(x(ti)) 



X+(*2) 



(6.226) 
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The reader should be aware that these expressions may require modification on the initial 
value surface to account for whatever operator ordering prescription is imposed upon the 
paraphernalia of surface gauge fixing. If both the time-ordered product of 0[x,7r] and 
the two states are invariant then relation (6.20) shows that (6.22a) and (6.226) give the 
matrix element and expectation value, respectively, of the reduced operator T^O[x r , ^ r ]j- 
If C[x,7r] is invariant before time-ordering then the necessary ordering corrections are 
the simple ones of the unconstrained theory, not the potentially complicated ones of the 
reduced theory. That is, we first convert the reduced operator 0[x r , 7r r ] into the analogous 
unconstrained operator, 0[x,7r], inside the matrix element or expectation value. Then we 
write it as the time-ordered product, t(o\x, tt] plus ordering corrections, and we use 
(6.22) to evaluate the matrix element or expectation value of the time-ordered product. 
Since the ordering corrections are also operators in the unconstrained theory we time order 
them and apply (6.22) again, repeating the procedure as often as is necessary to reduce 
everything to functional integrals. 

If the states and the operator in (6.22) are invariant then these expressions are inde- 
pendent of the residual gauge condition G k = 0. However, even when the wavefunctions 
or the functional C[x,7r] are not manifestly invariant the right hand sides of (6.22a) and 
(6.226) still represent the matrix elements or expectation values of some invariant operator 
in the presence of some invariant states. That this must be so follows from the fact that 
the gauge has been completely fixed and any quantity becomes invariant when it is defined 
in a unique gauge. To find which invariants one employs the residual gauge condition to 
"invariantize" the wavefunctions and the operator [21,26]. 

It is to some extent pointless to expend effort in order to discover a manifestly invariant 
state or operator when the only way physical information can be extracted from these ob- 
jects is by taking gauge fixed inner products. The aesthetic advantage to using manifestly 
invariant wavefunctions and operators is that then matrix elements and expectation values 
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are independent of the choice of gauge. If we understand physics in a particular gauge - 
as we often do — then this is not much of an advantage. The practical advantage is that 
manifest invariance allows us, through relation (6.20), to compute the matrix element or 
expectation value of a reduced operator 0[x r ,7r r ] using the same matrix element or ex- 
pectation value of the unconstrained operator 0[x, tt]. We can therefore avoid the need 
to ever construct the reduced Hamiltonian which, as we have seen, can be a formidable 
task. Of course we can always invariantize non-invariant states and operators, but then 
the non-locality and complicated field dependence this typically entails makes it difficult to 
compute the time ordering corrections. However, we emphasize that the process is simple 
enough to carry out perturbatively — witness section 4 — and the very fact that only 
operator ordering corrections are needed to relate C[x,7r] to 0\x r ,ix r \ inside gauge fixed 
inner products shows that we are well and truly free of the formal paralysis described as 
the paradox of second coordinatization. 

We turn now to a survey of how the general formalism manifests itself in the three 
dynamical systems we have studied. For the coupled oscillator model of section 2 the 
x a (tys are q\{t) and q2(t) while the 7r a (t)'s are p\{t) and P2(t). The evolution functions 
X a (x, 7r, and Tl a (x, tt, tj can be read off from (2.34). The unconstrained Hamiltonian 
is: 



Because it is quadratic in the momenta we can convert the canonical functional integrals 
of (6.22) into the more familiar, configuration space functional integrals. 

Recall that since the coupled oscillator model is not really a gauge theory we imposed 
P2 = as an ersatz "constraint." We might therefore identify the single constraint by 
inverting (2.34) to give p2 in terms of the time evolved variables: 




(6.23) 




(6.24) 
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The most general invariant operator is a function of q\, p\ and p 2 . We can of course express 
these initial value operators in terms of the evolved operators. The result for p 2 is just 
(6.24); the results for q\ and p\ are: 

[?l(*)-92(*) 

/3 \ 1 
ajt +_ 

V2 / ma 
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,'3 \ 1 
cos(-^j+- L 



cos ( -ut 



+ 



3muj 
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(6.25a) 
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(6.256) 



"Invariant" states are those whose wavefunctions are independent of q 2 in the t = basis 
of position eigenkets. Since invariant operators are independent of q 2 it is clear why any 
P2's can also be neglected when acting on an invariant state. 

Since the residual gauge condition is G = q 2 the associated Faddeev-Popov matrix is: 



M = -1 



(6.26) 



The single canonical pair (q,p) is of course just (qi,pi); the pure gauge variable is g k — ► 
q 2 , and the constrained variable is C£ — ► p 2 . Note that the unconstrained operators %(£) 
and Pi(t) are not invariant because they depend upon q 2 . This means that products of 
them do not give products of the analogous reduced operators in matrix elements and 
expectation values. For example, it is straightforward to show that in this gauge: 

^ sin (-ws) - sin ^ws) ^ cos (J^') - ^ cos (7^/) ^ inv ; t ^[ nv ; t' 



G 
(6.27) 



The terms on the extreme right represent unphysical and gauge dependent information that 
is included in the expectation value or matrix element of the unconstrained, non-invariant 
operator. 

For scalar QED in temporal gauge the x a (t)'s are the fields <J)(t, x), (j)*(t, x) and Ai(t, x); 
the 7r a (t)'s are ir(t,x), ir*(t,x) and Ei(t,x). Owing to the interaction, there is no closed 
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form expression for the evolution functions X a {x, tt, tj and Ii a ^x, tt, tj . The uncon- 
strained Hamiltonian is some ordering of (3.6). Since it is quadratic we can convert the 
canonical functional formalism into the configuration space version. 

The constraint — of which there is one for each space point x — is some ordering of 
(3.8). A typical invariant operator is: 



Cinv = <t>{t, x) exp 



te 



d 3 x'G(x;x'^ diAi(t,x) 



x Ej(u, y)ir(v,z) exp 



-ie J d 3 z G (^z; z d^A^{v,z) 



(6.28) 



where G^x;x is the Green's function first used in (3.13). The exponential term associ- 
ated with each of the charged fields creates the associated longitudinal Coulomb field. A 
typical invariant wavefunctional is that of the free vacuum which, up to a normalization 
factor, is: 



mv 



oc exp 



d 3 x$*(x) v /Z V2$(x) 



<&(x) = 4>{x) exp 



-\jd\ Fl] {x)^L^F, 
ie J d 3 x' G^x; x '^j di Ai(x) 



(x) 



(6.29a) 



(6.296) 



The apparent coupling between the vector potential and the charged fields is a gauge 
fiction; it disappears when the instantaneous Coulomb gauge condition (3.11) is used. 
This is an example of why it is pointless to struggle to achieve manifest invariance in an 
expression which must ultimately rest against a gauge fixing delta functional. Note as well 
that (6.28) and (6.29) are only invariant with respect to the residual, time independent 
gauge symmetry. 



For instantaneous Coulomb gauge, (3.11), the Faddeev- Popov matrix is: 



M 



(x;yj 



•TT^TT"^ 3 (x - y) 

ox 1 oy l V / 
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(6.30) 



Note that it has a non-singular determinant on the appropriate function space. The most 
natural identification of canonical, gauge and constrained operators on the initial value 
surface is: 

q a ^{AJ(x)^(y)^(z)} (6.31a) 
p a ^{E?(x)My)S*(?)} (6-316) 



Cfc 



Af(x) 
Ef{x) 



(6.31c) 



(6.31d) 



Note that although these initial operators do commute canonically — on the appropriate 
function space — we saw in section 3 that their time evolved versions do not. If we assume 
that t > u > v > then the invariant operator (6.28) is also time ordered and we can 
write: 



vl>. • f 



TO 



*in V ;0) G = 



dAu(s,w) 4>(s,w) d(f>*(s,w) S Aq(s,w) 
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-ie I d 3 z' 'G[z;z')d k A k (v,z) 



(6.32) 



x exp 



t 



i i ds I d 3 wC 



d t A e (0, w)] * inv [A(0, w), 0(0, w), <f>*(0, w)] 



where expression (3.1) defines the Lagrangian C. Note that we have integrated out the 
momenta and that we have introduced Aq as an integral over the volume gauge fixing delta 
functional. Since the operator (6.28) is invariant, as are the states, it follows that this 
functional integral represents the matrix element of the corresponding reduced operator. 
The interested reader can worked out examples from standard QED in ref. [26]. 

For quantum general relativity at fixed lapse and shift the x a (tys are the 3- metrics, 
7jj(£, x); the n a (tys are their conjugate momenta, tt 1 ^ (t, x). It is no more possible than for 
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scalar QED to give explicit forms for the evolution functions, X a (x, tt, tj and n a (x, tt, tj . 
The Hamiltonian is some operator ordering of the classical result: 



H 



7,7T 



(t) = J d 3 x 7, tt (£, x) Hn 7, tt (t, x) (6.33) 



Note that the lapse and shift may depend upon time and also upon the dynamical variables; 
in fact dependence upon the dynamical variables is necessary classically in order to have 
a chance of avoiding the evolution of coordinate singularities. Although no one has ever 
exhibited a gauge which is classically free of coordinate singularities its existence seems 
obvious if a sufficient amount of field dependence and non-locality is permitted in the lapse 
and shift. In any case we shall assume that such a gauge exists. 

Note that if the lapse and shift do depend upon the 7r^'s then the canonical action may 
not be quadratic in the momenta and we would not be able to reach the usual configuration 
space functional formalism. Although inconvenient, there is no inconsistency in this. It 
could happen as well in scalar QED if we had allowed Aq to depend upon n or n*. Note 
also that even if the lapse and shift depend only upon 7^, integrating out the momenta 
will give a field dependent measure factor. If the lapse depends non-locally upon 7^ then 
this measure factor can be non-local, and it may therefore give non-zero contributions even 
in dimensional regularization. Again, this is not specific to gravity; we could have made it 
happen as well for scalar QED by permitting Aq to depend non-locally upon the fields. 

There are four constraints for each space point x. Without regard to operator ordering 
they are:* 

2 

H = ^— (~f lk ljt - ^7ij 1m) ^ 3 k M - ^ (R - 2A) ^ (6.34a) 



n % = -2 ll3 ^ k . k (6.346) 



* It is futile to pay much attention to ordering so long as the problem with renormalization 
precludes being able to take the unregulated limit of inner products [27]. 
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It is reasonably straightforward to construct a basis of operators which transform as scalar 
densities under an arbitrary diffeomorphism. The method is to evaluate Riemann tensors 
at the ends of geodesic segments emanating from a central origin, with the indices of each 
Riemann tensor defined in the local inertial frame field of the origin as obtained by parallel 
transport along the connecting geodesic [28]. If the coordinate time generated by (6.33) 
can be extended infinitely then full diffeomorphism invariance results when the origin is 
integrated over R x M 3 . Since such operators are invariant under all diffeomorphisms they 
are necessarily invariant under those of the residual symmetry group which preserve the 
fixed lapse and shift. 

No one has obtained closed form expressions for normalizable invariant states, but one 
can of course construct them order by order in perturbation theory. It is the fact that the 
asymptotic functional formalism does this automatically which explains the observation in 
[26] and [28] that all of Einstein's equations are obeyed in regulated perturbation theory. 
Note also that there is a considerable phenomenological advantage to building invariant 
states perturbatively: one knows what they mean this way. If given a solution of the 
constraints it is difficult to tell what level of excitation it represents since on a closed 
spatial manifold all states are degenerate with zero energy. However, if this solution 
is perturbatively related to the invariantized free vacuum then one can assume - - at 
least as a first guess — that it represents empty space. Similar comments apply to the 
wavefunctionals built up from other Fock space states. 

In section 4 we reduced the theory in two steps. In the first step we enforced the con- 
straints that fixed and pi and surface gauge fixed their conjugate variables. Explicitly: 

The Cfc's are: 



h % = + C(k 1 ) 



(6.35a) 



p i = + O(K 1 ) 



(6.356) 
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while the G^s are: 



p* = 



hi = 



(6.36a) 
(6.366) 



The Faddeev-Popov determinant arising from this first step of reduction is just a C-number 
to lowest order. This is not the case for the second step. In the second step we enforced 
the constraint on p tr and surface gauge fixed h tr using: 



V 



1 =6Eq + 0(k 1 ) (6.37a) 

h tr = (6.376) 

We see now that the Faddeev-Popov matrix for this, the second step is 2p tr . The choice 
of variables on the initial value surface goes as follows: 

q a — >h tt (6.38a) 



Pa 



fit 



(6.386) 

f^{fi,h\h tr } (6.38c) 

c fc — > {S*,^} (6.38a 1 ) 

We will see next how the inner product breaks up in two parts. One part for negative p tr 
and one for positive, representing and expanding and a contracting space respectively. 

Let us look at the inner product defined in (6.18): 



^25*2 



4>l\tl ) = (^25*2 
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G[x 

tr 



abs < 2 p 
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(6.39) 



Where we used ~ because we are disregarding overall multiplicative factors and we are 
working to lowest order in k. As mentioned in section 4 each wavefunction is divided in 
two parts depending on the action of p tr upon them: 
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(6.40) 



The inner product of (6.39) can then be written as: 



^2^2 



^r,h ) = (^2;*2 



e(V r ) ^(V r ) p <r e(V r ) - e(V r ) P* r ^(/i <r ) e(V r ) 



(6.41) 



Where we have chosen a Hermitian ordering. Using (6.40) equation (6.41) becomes: 



^25*2 



+ (^2 ; £ 2 



(6.42) 



We see that, as previously promised, the inner product breaks up into two parts. Both 
of these parts are present quantum mechanically; however, classically either ip + = or 

i>- = o. 



7. Conclusions 

The context of our analysis is that of a gauge theory in which the ability to perform 
local, time dependent transformations has been fixed but there is still a residual gauge 
symmetry characterized by how it acts on the initial value surface. Examples of such 
theories include scalar electrodynamics in temporal gauge, which was studied in section 
3, and general relativity with fixed lapse and shift, which was studied in section 4. The 
initial value problem has a unique solution in this setting — which it does not in the 
fully invariant theory — because residual symmetry transformations cannot change the 
time evolved dynamical variables without also changing their initial values. We assume 
that a canonical formalism describes this initial value problem; we assume as well that 
the residual gauge symmetry is generated in the usual way by functions of the dynamical 
variables which are constrained to vanish as a consequence of the Euler-Lagrange equations 
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of the (volume) gauge fixed variables. This is what we mean by the "unconstrained theory," 
and its generic dynamical variables are the x a (t)'s and 7rg(t)'s of section 6. 

The generic reduced theory is obtained by fixing the residual gauge freedom, imposing 
the constraints, and identifying a subset, <jV(t)j, of the original dynamical variables 
which gives a complete and minimal representation of physics. What this means is that the 
constraints, the residual gauge conditions and the values of the 's at any time t uniquely 
determine the a; a (t)'s and the 7rg(£)'s. The reduced dynamical variables inherit their time 
evolution and their equal time bracket (or commutation) algebra from the unconstrained 
theory. This is true even when the unconstrained Hamiltonian vanishes after reduction; it 
is even true if no Hamiltonian exists which generates the time evolution of the w z (t)'s. 

Section 2 described a standard construction of the last century for identifying canonical 
variables q a (t) and p^(t) in the reduced theory, and for finding the non-zero Hamiltonian 
which generates their time evolution. The identification of a reduced canonical formalism is 
not unique. On the classical level we can vary it by performing canonical transformations; 
with suitable attention to operator ordering there is a similar class of transformations on 
the quantum level. What keeps physics unchanged is the fact that measurement theory is 
based upon the x a (t)'s and the 71^3 (t)'s, considered as functions of the reduced canonical 
variables. As our identification of the reduced canonical variables changes, the functional 
dependence of the unreduced canonical variables changes so as to keep the x a (t)'s and 
7Tg(t)'s the same. 

An immediate consequence of the multiplicity of reduced canonical formalisms is that 
the associated Hamiltonians have no physical significance independent of whatever meaning 
attaches to the reduced canonical variables whose evolution they generate. Although this 
may seem strange, quantum field theorists are familiar with the phenomenon through 
the "interaction representation." This is a field redefinition in which any perturbatively 
well defined interacting quantum field theory is transformed into the corresponding free 
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quantum field theory. We do not conclude that all perturbatively well defined quantum 
field theories are free because we insist that physics be inferred from the original fields. 

The construction of a reduced canonical formalism is irrelevant to most issues in clas- 
sical physics. This is because we infer physics from the x a (£)'s and 7Tg(t)'s, and we may 
as well solve for these variables directly in the unconstrained theory, starting from initial 
value data which obey the constraints. The principle motivation for erecting a reduced 
canonical formalism is as a prelude to applying canonical quantization. This is especially 
relevant to systems, such as gravity on a spatially closed manifold, for which reduction 
causes the unconstrained Hamiltonian to vanish. As has been noted, we did not invent the 
construction discussed in section 2; this was done by the classical physicists of the previous 
century [5,6,10-15]. Our contribution is rather to propose that quantum gravity should 
be defined by canonically quantizing a reduced canonical formulation of whatever is the 
correct theory of gravity. 

We constructed explicit reduced canonical formalisms for the coupled oscillator model 
of section 2, and for scalar electrodynamics in section 3. Although no complete construc- 
tion seems possible for general relativity we showed in section 4 how to do it perturbatively 
around a flat background on T 3 x R. We also gave explicit constructions for a handful of 
minisuperspace truncations in section 5. Our inability to give an explicit, non-perturbative 
construction for the complete theory of general relativity does not pose a practical barrier 
- even for the study of non-perturbative phenomena — because of the relation we were 
able to obtain in section 6 between the relatively simple quantum mechanics of the uncon- 
strained theory and that of the reduced theory. Equation (6.20) asserts that the matrix 
elements or expectation values of invariant functionals of the reduced operators are equal 
to the matrix elements or expectation values of the same functionals of the unconstrained 
operators in the presence of invariant states. Explicit functional integral representations 
exist for unconstrained matrix elements and expectation values. These representations are 
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as simple to evaluate as it seems possible for anything to be in an interacting quantum 
field theory, and they have at least the potential for extension beyond perturbation the- 
ory. Thus it is not really necessary to explicitly construct a reduced canonical formalism; 
it suffices to know that such a construction exists and that we can study it using the far 
simpler formalism of the unconstrained theory. 

It is worth reviewing how the formulation we propose for canonical quantum gravity 
avoids the four correspondence paradoxes mentioned in section 1. Our resolution to the 
paradox of second coordinatization is that the lapse and the shift determine what is meant 
by time evolution in quantum gravity the same way they do for classical gravity.* It is 
neither necessary, nor even particularly desirable, to try doing the job twice by identifying 
some other variable as "time" and then attempting to interpret the Wheeler-DeWitt con- 
straint as a Schrodinger evolution equation. Our method works because the Heisenberg 
field operators depend upon the time fixed by the lapse and the shift, whether or not we 
restrict the initial value data by fixing the residual gauge and imposing the constraints. 
The conventional method will only work if one of the Heisenberg operators is an invertible 
function of time, and the method will only produce tractable results if this time depen- 



A tangential point concerns the order of gauge fixing. Our volume gauge condition de- 
termines the lapse and shift as functionals of Jij(t, x), n l] (t, x), and possibly also of time 
and space; we then impose a surface gauge condition upon the 7r/'s and 7r^'s, and use the 
constraints to determine the evolution of some components of the 7jj's and 7r^'s. Many 
researchers [3,7,8,9] prefer to impose a volume gauge condition upon 7^ (t, x) and 7r z -?(£, x); 
they then use the constraint equations, with some surface gauge condition, to solve for the 
lapse and shift. An example of our method in scalar electrodynamics is fixing Ao(tx) = 
as the volume gauge condition, and then using the constraints and the surface condition, 
diAi(0,x) = 0, to determine the longitudinal field components. An example of the other 
procedure would be fixing c^A^t, x) = as the volume gauge, and then using the constraint 
equation plus the freedom to perform time dependent, harmonic gauge transformations, 
to determine Ao(t, x). Our method gives the best chance of defining a successful evolution 
in gravity because it allows one to adjust the rate of evolution in response to what the 
fields are doing. In this way we can avoid coordinate singularities, and also — by merely 
slowing down before they form — the evolution of physical singularities. Of course if the 
other method gives a successful evolution it can always be transcribed into our form by 
merely regarding the derived solutions for the lapse and shift as volume gauge conditions. 
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dence is relatively simple. Gravity is a formidably complex dynamical system and it may 
somewhere harbor such an operator; it is certain though that no one has found it yet, 
despite years of search. 

A key point of our proposal is that one should infer physics the same way in the 
quantum theory as in the classical theory: by studying the reduced variables, x^(t) and 
7Tg(£). The fact that these variables are not manifestly invariant before gauge fixing is 
irrelevant. Any quantity can be made gauge invariant by defining it in a particular gauge 
and this is just what fixing the residual gauge freedom does. No classical relativist would 
demur from modeling planetary orbits using the geodesies of the Schwarzschild metric just 
because the latter was reported in Schwarzschild coordinates. Our formalism permits the 
same indifference on the quantum level. Correspondence in the limit that h vanishes is 
manifest through relation (6.20) and the functional expressions (6.22) for unconstrained 
matrix elements and expectation values. 

Our resolution to the paradox of dynamics is that proper correspondence should pertain 
in the limit that G vanishes if that the reduced canonical variables are taken to include the 
canonical variables of the pure matter theory. In this case the Hamiltonian of the reduced 
canonical formalism for gravity plus matter will go over into that of free gravitation around 
some background, plus the pure matter Hamiltonian for that background. The same 
correspondence will obviously apply for inner products. We have no general proof that a 
suitable choice exists for the reduced canonical variables but, if it does, then the result 
follows automatically from the correspondence limit of the evolution equations.* We feel 
it is plausible that such a choice exists, and we have given an explicit example in the 
minisuperspace truncations of section 5. 

Our resolution to the paradox of topology is that a correspondence limit should exist 
between the non-zero Hamiltonians of infinite, spatially open manifolds [3,4] and the re- 
* That the field equations obey the correct correspondence limit is not in doubt. 
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duced Hamiltonians of spatially closed manifolds. This correspondence will exist when the 
reduced canonical variables are chosen to include the canonical variables of the open space 
theory, and when the coordinate volume of the closed space is taken to infinity in such a 
way that the initial value data obey the asymptotic condition of the open space. The idea 
is that in such a situation causality prevents the closed topology from reaching around the 
universe to influence the localized initial disturbances. As with the paradox of dynamics, 
the desired correspondence follows immediately from the field equations provided that a 
suitable choice can be made for the reduced canonical variables. We have no general proof 
that such a choice exists but we did show in section 4 that it does for gravity on the 
manifold T 3 x R. 

We comment that if the aforementioned topological correspondence limit holds gener- 
ally then the 2 + 1 dimensional constructions of Moncrief [7] , Hosoya and Nakao [8] and 
of Car lip [9] would be special cases of the general formalism described in section 2. The 
same would be true of the 3 + 1 dimensional constructions of A.D.M.[3], and of Deser and 
Abbott [4]. Indeed, the method of section 2 seems to provide the long sought unifying 
principle needed to define energy on a space of arbitrary topology. 

We do not avoid the paradox of stability by appealing to the fact that reduced Hamil- 
tonians have non-trivial spectra. As has been noted, what the reduced Hamiltonian is 
depends upon how we identify the reduced canonical variables, so its spectrum for an ar- 
bitrary identification has no intrinsic significance. Nor are reduced Hamiltonians typically 
conserved. They do not therefore fill the usual role of constraining the way in which a 
theory can excite its various degrees of freedom. In fact for gravity on a spatially closed 
manifold the only generally conserved energy functional is the unconstrained Hamiltonian 
which vanishes upon reduction. This means that all states really are degenerate and hence 
that the universe is liable to evaporate into pairs in the manner of section 1. 

That we are here is a consequence — assuming that our spatial manifold is closed - 
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of causality and of the weakness of the gravitational interaction. The H = constraint 
is not met, as is sometimes supposed, by a vast reduction in the number of possible 
states compared with gravity on an open space. The perturbative analysis of section 4 
shows that there are at least as many positive energy graviton modes on T 3 x R as on 
i? 3 x R; in fact the H = constraint is enforced by the global negative energy mode, 
p tr . We conjecture that this is the case generally; that is, positive energy modes can 
only be excited by corresponding excitations of a global negative energy mode. But a 
global mode by definition pervades the spatial manifold, so exciting it requires a similarly 
extensive process. On large manifolds causality imposes a formidable barrier to such 
excitation. There is an additional barrier in the fact that the global mode can only be 
excited gravitationally. Since gravitational interactions are typically very weak in our 
current universe they must proceed slowly. Note that neither barrier would apply to a 
strongly gravitating system of small physical volume. One consequence of our work is the 
prediction that such systems ought to be unstable. 

The barrier causality imposes against instability becomes absolute in the limit that the 
coordinate volume goes to infinity while the initial value data are only locally disturbed 
from a vacuum solution. In this limit the global mode decouples, and both its (negative) 
energy and the (positive) energy of the local modes become separately conserved. This is 
how we can approach the conserved Hamiltonians [3,4] of spatially open manifolds. The 
gravitational barrier becomes similarly infinite in the limit that G vanishes. In this limit 
the negative mode again decouples — along with all the other gravitational modes — and 
the energies of matter and of each gravitational mode become separately conserved. This 
is how we can approach the conserved Hamiltonians of pure matter theories. 

We should comment as well on a venerable argument sometimes used to deny the 
possibility of non-zero Hamiltonians for gravity on closed spatial manifolds. The argument 
begins with the observation that any such Hamiltonian would have to be the integral of 
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a local function of the metric and its first derivative. If it is to have physical significance 
such a Hamiltonian must also be an invariant, but there are no invariant functions of the 
metric and its first derivative. Therefore the Hamiltonian of any theory with dynamical 
gravity must be zero when the spatial manifold is closed. We evade this argument by 
the simple device of permitting the Hamiltonian to be non-invariant. It generates time 
evolution in a particular coordinate system — the one provided by the fixed lapse and 
shift — so it can and should be a gauge dependent object. The physical energy, on the 
other hand, must really be an invariant, and we have seen that it is zero. There is no 
contradiction between these two facts because the non-zero Hamiltonian that generates 
time evolution for a reduced dynamical system need not be the physical energy, neither 
must it be necessarily conserved. We cobbled together a simple example of this with the 
coupled oscillator system of section 2, and the phenomenon is generic in gravity. 

We close with an admonition to those who would prefer to continue searching for a 
conventional solution to the "problem of time" — i.e., the identification of some component 
of the metric as physical time and the interpretation of the Wheeler-DeWitt constraint as 
a Schrodinger evolution equation. Though our frank admission that time is a gauge choice 
and our frequent appeal to perturbation theory in implementing this gauge choice may 
seem ugly, the construction of section 2 is not intrinsically tied to perturbation theory - 
witness the functional representations of section 6 — and it has the great advantage that 
one can tailor the coordinate system to the state and the operator under study. There is no 
reason why a "magic bullet" should even exist for identifying physical time, much less that 
it should be a simple functional of the metric. Nor is it in accord with our experience in 
other branches of physics to assume that a single calibration should govern all conceivable 
measurements. We do not insist, for example, that condensed matter experimentalists use 
the same thermometers in the microkelvin regime as in the heart of an exploding, high 
field magnet. Nor do we deny the validity or utility of measuring temperature just because 
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different techniques are employed in different regimes. What we insist upon is rather that 
any two methods should agree in those environments for which they both apply. There 
is no reason to impose a more stringent standard in the vastly more varied and extreme 
environments imaginable to a theorist in quantum gravity. 
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